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Chapter 1 

Introduction 



Elastic light scattering (LS) has been extensively used to study samples showing 
a non uniform refraction index on lengthscales from a fraction of a microme- 
ter to a fraction of a millimeter. Basically, a wide laser beam is sent through 
the sample, and the light scattered at any angle is measured by a detector in 
the far field. Many phenomena have been studied with this technique; among 
them, thermodynamical fluctuations of concentration in solutions, of tempera- 
ture, of pressure, convcctive instabilities and turbulence, colloids and colloidal 
gels, systems showing critical opalescence. 

The measure of the intensity of the beams scattered by a target is a well 
known way for analysing physical properties of a sample, and is used in many 
fields of physics: examples are X ray diffraction in crystallographic analysis and 
particle scattering by nuclei in the Rutherford experiment. In these examples, 
as in light scattering, the intensity of the scattered beams gives informations 
on the square modulus of the Fourier transform of a quantity of interest of the 
sample, the so called power spectrum, evaluated in the transferred wavevector. 
In the case of light scattering, the quantity of interest is the refraction index. 
Since refraction index shows variations due to concentration, temperature and 
pressure fluctuations, all these quantities can be investigated using light scat- 
tering. 

In principle, measuring the intensity of the light scattered at an arbitrary 
small angle allows to obtain informations on arbitrary long wavelength Fourier 
components of the sample. In facts, the measurement of the intensity of the 
scattered light becomes more and more difficult as the scattering angle becomes 
small, mainly due to the stray light, that is the light scattered by the imperfec- 
tions of the optical system, which is mainly scattered at small angles. Thus light 
scattering cannot give informations on the features of a sample, if their associ- 
ated wavelenght is longer than a given value. On the other hand, image forming 
techniques, such as Schlieren, dark field, phase contrast microscopy, have no 
limitations on the size of the features they can observe. The limitation of many 
image forming techniques is the difficulty to quantitatively relate the observed 
images to the physical properties, mainly when dealing with three dimensional 
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samples. 

In the present work, we describe three new techniques, which aUow to mea- 
sure the scattering intensities, overcoming the difhcuhics associated to small an- 
gles. The first of these techniques, the hOmoyne Near Field Speckles (ONFS), 
has been presented very recentely [T]; in the present work, we show new re- 
sults, obtained with a slightly improved optical setup [5]. The other techniques, 
the hEterodyne Near Field Speckles (ENFS) and the Schlicren-like Near Field 
Speckles (SNFS) are improvements based on that. The first has been recently 
patented by us [21 U; the second is presented here for the first time. Moreover, 
in Chapt. [31 we present for the first time a mathematical derivation of the 
working formulas for the three techniques. 

Basically, the experimental setup consists in a wide laser beam passing 
through the sample; a lens forms an image of a plane at a given distance from 
the cell on a CCD sensor. The image, in the near field, shows speckles, since 
it is formed by the stochastical interference of the light coming from a random 
sample: the electric field has a gaussian probability distribution. We will show 
that, under suitable conditions, the correlation function of such a field closely 
mirrors the correlation function of the investigated sample; moreover, in general, 
from the correlation function of the speckle field we can calculate the scattered 
intensity. 

The lens that forms the image on the CCD focuses the transmitted beam 
around a given point. In ONFS, a beam stop is placed in that point, in order to 
dispose of the transmitted beam; in SNFS a blade stops half transmitted beam, 
along with one half of the scattered light, like in Schlieren technique; in ENFS no 
opaque element is introduced in the optical system. In ONFS, the CCD sees the 
speckles given by the interference of the scattered beams with themselves: ONFS 
is an homodyne technique . In ENFS and SNFS, the speckle field is heterodyned 
with the much more intense transmitted beam, that acts as a reference beam: 
the measured intensity is linear in the speckle electric field. We acquire a set of 
images, from one of the three techniques, by using a CCD camera, connected 
to a frame grabber; the images are then elaborated by a PC, to obtain I{q), 
the scattered intensity as a function of the transferred wavevector, the same 
information obtained by LS. For each technique, we developed algorithms which 
allow to evaluate the scattered intensities, including the corrections for the stray 
light and for the shot and read noise of the CCD camera. The algorithms are 
described in Chapt. [5]and[6l 

We used ONFS and ENFS to measure the scattered intensity of some col- 
loids, and wc compared the results with those made by a state-of-the-art clas- 
sical Small Angle Light Scattering (SALS) device. The agreement is very good, 
notwithstanding a much simpler and stable layout. In Chapt. [7] we compare 
ONFS, ENFS and SALS measurements, and discuss the main sources of errors. 
Scattering intensities measured with ENFS show a better quality; we used this 
technique to evaluate the diameter distribution of some known colloids, by using 
an inversion algorithm based on Mie theory. The results are presented in Chapt. 
[H this shows that ENFS is a simple and powerful alternative to SALS, suited 
for industrial applications like particle sizing. Moreover, SNFS has been used 
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to evaluate the power spectrum of non-equilibrium fluctuations in a free diffu- 
sion process, thus showing that such techniques have interesting applications in 
fundamental physics; results arc shown in Chapt. [51 
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Chapter 2 



Qualitative description of 
the technique. 

The intensity of the hght scattered from a spatiaUy disordered sample has a 
speckled appearance, the speckles being generated by the random interference 
of the scattered elementary spherical waves. While the study of the one point 
intensity time correlations has proven very useful, and it has generated the 
technique of Intensity Fluctuation Spectroscopy (IFS) [5], the measurement of 
the two point, equal time, intensity space correlation function, that is the size 
and the shape of the speckles, does not provide any useful information. Indeed 
the Van Cittert and Zernike theorem states that the far field space correlation 
function depends only on the intensity distribution of the scattering volume, 
and in no way depends on the physical properties of the sample. 

In this chapter we will present qualitative elements showing that for fluctu- 
ations the size of the wavelength of light or larger, in the near field we obtain a 
speckle field, that is, a gaussian field; moreover its statistics is directly related 
to the scattered intensity distribution. We will derive the working formulas 
for three tecniques, hOmodyne Near Field Speckles (ONFS), hEterodyne NFS 
(ENFS) and Schheren-like NFS (SNFS); analogies with the IFS wiU be pointed 
out. Advantages with respect to the more conventional Small Angle Light Scat- 
tering (SALS) technique will be discussed. 

First of all, we will describe ONFS setup; many considerations hold also for 
ENFS and SNFS. The experimental set-up is very imorthodox, with respect to 
a conventional SALS device. It consists of a wide laser beam and of a Charge 
Coupled Device (CCD) detector positioned so to be flooded with light coming 
from any scattering direction the system can scatter at. 

The Van Cittert and Zernike theorem states that the field correlation func- 
tion is [B]: 



CEiAx, Ay) = (E (x, y) E* {x + Ax,y + Ay)) 



I{^,ri)exp 



{^Ax + ijAy) 



d£,dr]{2.l) 



where E {x, y) is the field in the observation plane x — y, A is the wavelength 
and / rj) is the actual intensity distribution of the source in the plane ^—77 at a 
distance z from the observation plane. The theorem holds for sources consisting 
of point emitters, like atoms. The intensity correlation function Ci (Ax, Ay) = 
(/ {x, y) I {x + Ax, y + Ay)) is then derived by applying the so called Siegert 
relation [7]: 



Ci {Ax, Ay) = {if + \Ce {Ax, Ay) p (2.2) 

Equations (|2.ip and (|2.2p specify that the intensity correlation function is 
related to the space Fourier transform of the source. In practice, this implies 
that a source of size D will generate speckles of size -^z on a screen positioned 
at a distance z . 

We will start introducing simple euristic arguments and crude evaluations 
for the near field speckles of the scattered light. Let us consider the case of 
a large beam diameter D, impinging onto a sample of particles of diameter d 
larger than the wavelength of light: see Figure [2?]T a). Most of the power will 
be scattered in a forward lobe of angular width « ^. Let us consider a small 
area S, for example a multi-element sensor array, in the immediate vicinity of 
the scattering volume: see Figure [^?lT b). Let us assume that we can ignore 
the transmitted beam: we will take care of this problem later on. Although the 
sample is illuminated over the entire surface of diameter D, the light falling onto 
the sensing area will come only from a smaller area of diameter D* . In fact the 
brightness of the scattering volume will change as a function of the observation 
angle in a way that mirrors the scattered intensity distribution. Consequently, 
for the sensing area, the source region from which light is drawn is a circle with 
a diameter D* = ^z, z being the distance of the sensing area from the scattering 
surface; source regions outside do not contribute appreciably. We say that the 
near field condition is met li D* <C -D. One can then immediately estimate the 
size of the speckles dsp = -^z fa d, a remarkable result in many respects! The 
speckles have the size of the particle diameter, and this value does not depend 
on the distance z from the sample, provided the near field condition D* <C D. 
This has to be compared with far field speckles, whose diameter scales linearly 
with the distance from the source. Also notice that the actual sample thickness 
does not matter, provided that the near field condition is met, and that the 
speckle size does not depend on the light wavelength, an unexpected feature for 
an interference pattern. 

Notice that all the above applies under conditions that arc more stringent 
than the usual "near field" condition [8] for a source of size D , namely -pi ^ 1 . 
In the present case the condition is D* <^ D which implies -53 ^ 1- 

To put things in a more quantitative way, we will determine the near field 
intensity correlation by first re-writing the Van Cittert and Zernike theorem in 
a more appropriate form. We notice that Eq. ()2.ip may be rewritten in the 
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Figure 2.1: (a) Small angle scattering. A beam of diameter D impinges onto a 
sample composed of particles of diameter d. Any zone within D will scatter light 
into a lobe of angular width Q = \/d (the length of arrows indicates scattered 
intensity), (b) Same sample, as in part (a). A sensor S close enough to the 
sample will draw light only from a zone of radius D* < D. Regions outside, 
even if illuminated by the main beam, do not feed light to S. Notice that again 
e = X/d. 
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following way: 



Ce{i^^ j Hq)e'^'Aq, (2.3) 

where r = (Ax, Ay), and g is a vector whose components are = and 
qy = j^T], which equals the scattering wave vector for small scattering angles. 

Equation (|2.3p is only a different way of writing Eq. (|2.ip . and I {q) is 
the intensity distribution of the source as seen from the observation plane as 
a function of the scaled angles {2it / X) / z) , and (27r/A) (77/0). As discussed 
in the introductory remarks, in the very near field / {q) equals the scattered 
intensity distribution, which is proportional to the Fourier transform of the 
sample density correlation function g (r) = {51 (r) 51 (0)), where 51 is the local 
fluctuation of the particle number density, integrated over the light path. Then, 
from Eq. it follows that: 



Cjir)^{lf[l + \girT\- (2.4) 

We would like to point out that Eq. (|2.4p closely duplicates the well known 
relation that holds for the IFS < I{0)I{t) >=< I [1 + \g(t)\'^] , where g{t) 
is the time correlation fmiction (see for example [S]). It should be noted that 
for some scatterers the Raylcigh Gans approximation is invalid, for example for 
larger spheres where the Mie theory applies, and therefore the pair correlation 
function g(r) cannot be extracted from the scattered light. It remains true how- 
ever that even in those cases the correlation method permits the determination 
of the scattered intensity distribution / (g), by using the relation: 



|2 



Ci{r) = {lf 



I (q) e'-J-^-dg 



(2.5) 



To determine the spatial intensity correlation of Eq. (|2.4p . one must first 
obtain experimentally the instantaneous intensity distribution of the near field 
scattered light. In order to evaluate the intensity correlation function with rea- 
sonable statistical accuracy it is also imperative to gather intensity distributions 
over a substantial number of points. To this end a CCD is ideal, the number 
of pixel being larger than 10''. As we shall see, it actually turns out that one 
frame is enough for a fair acquisition of the correlation function. 

In a previous work [T] , some measurements have been performed on a scat- 
tering model, an opaque metallic screens with pinholes of 140 and 300 microns 
chemically etched in random positions. The surface fraction occupied by the 
pinholes was around 10% and 20% respectively. Experimentally this greatly 
simplifies the problem, since the scattered field is stationary and also there is 
no transmitted beam. We call this configuration hOmodyne Near Field Speck- 
les, since the signal is given by the interference of different scattered beams. 
Being a two dimensional sample, the scattered intensity was simply related to 
the correlation function of the transparency function T {x,y) with T = 1 inside 
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the pinholes and zero outside [B]. A Helium Neon parallel beam with diameter 
(jT points) D = 15mm was sent onto the samples, and the speckle field was 
recorded with a CCD at various distances z ~ 50cm, z — 75cm and z ~ 100cm 
0. The corresponding values for D* ranged from 1mm to 4.3mm so that the very 
near field condition was always met. The rather large dimension of the pinholes 
was chosen so that the speckles were appreciably larger than the CCD pixel size 
(typically 9/im). For each type of pinholes, the measurements performed at the 
three distances showed minute differences. The results are shown in Fig \2.2\ 
where the data are compared with the correlation functions of digitised images 
of the set of pinholes on the metallic screen H. Since in this case the sample is 
two-dimensional, g {r) is the correlation function of T {x). The width and shape 
of the main peak are fairly well reproduced, in spite of the limited number of 
frames used (four frames on statistically equivalent samples for each type of 
pinholes). 

While the data obtained with the screens prove that near field speckles do 
mirror the properties of the scatterers, we feel that to assess the desirability 
of the technique for realistic applications (for example in colloid physics) mea- 
surements had to be taken with particle solutions down in the micron range. 
In order to do this, three problems had to be solved. The speckles in the near 
field close to the cell have dimensions around one micron and therefore are too 
small for the available CCD pixel size. Also, one must dispose of the transmit- 
ted beam. Finally, the speckle intensity distribution must be frozen at a given 
instant. 

The first two problems have been solved with the simple optical arrangement 
shown in Fig. 12.31 A wide parallel beam is sent onto the sample, placed against 
a large aperture lens of focal length /. A wire is stretched in the focal plane to 
intercept the main beam. The CCD is placed a distance q away from the lens. 
The system magnifies the speckle size by a factor M = {q~ f)/ f > 1, and the 
scattering angles are decreased accordingly. It is illuminating to point out that 
the technique can be considered as a scaled down version of the classical Hanbury 
Brown and Twiss [lOj experiment where the star intensity distribution is mim- 
icked by the scattered light intensity patch in the focal plane, and the ground 
based intensity correlation arc the CCD intensity correlations. The unavoidable 
presence of the lens and its finite aperture introduces some complication with 
respect to the lensless arrangement used for the pinholes. In a previous work 
[1], the poor numerical aperture of the lens introduced a non-uniform trans- 
fer function, which had to be evaluated by measuring a known sample. In the 
present work, care has been taken in order to avoid such problems. 

When the scattered speckles are observed with the CCD in real time, one 



^For very large objects the near field is not really near, since one must let the difl^raction 
figures from various objects interfere. This leads to the additional condition D* = ^ 5, 
where S is the typical distance between scatterers. 

^Experimental correlation functions shown in figure arc calculated according to the follow- 
ing procedure. First, the Fast Fourier Transform of the intensity distribution is calculated. 
Then the result is squared. The auto-correlation is then obtained by anti-transforming the 
power spectrum. 
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Figure 2.2: Measured intensity auto-correlations as a funetion of displacement r 
for two sets of randomly positioned pinholes, of 140/xm and 300/im in diameter. 
For both the samples, measurements at three distances are reported, together 
with \g (r)l^, calculated from the digitized images of the two samples. 
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Sample 




Figure 2.3: Optical layout for ONFS. The main transmitted beam is blocked by 
a stop in the focal plane. Almost all the scattered light is sent to the CCD. 



notices quite vividly that the speckle size changes as the size of the scattcrers is 
changed. Also, for a given sample the speckles boil with the same time constant 
on the whole screen, the time constant getting larger for samples with larger 
diameter particles. With regard to the third problem mentioned above, these 
observations also indicate that even with a conventional CCD and a small power 
He-Ne laser there is no problem in getting instantaneous pattern distributions. 
Indeed even for the smallest particles that can be studied with present experi- 
mental set-up, with diameters down to lfj.ni, and assuming diffusive motion, the 
shortest time constant associated to the smallest scattering wavevcctor yields 
Trnin = 0.125s, a time long compared with the shortest frame exposure available 
with standard frame grabbers, typically l/16000s. 

Let us compare the Near Field Speckles technique with the more traditional 
Small Angle Light Scattering. The essential feature of a scattering layout [TTlfT^ 
is that the light scattered at a given angle hits the sensors along a circle of given 
diameter around the optical axis. We believe that the correlation method of 
NFS offers some distinct advantages over the scattering technique. First, there 
is no need for accurate positioning of the CCD, that can be rather casually 
placed at a distance z from the focal plane (see Fig. 12. 3p . At variance, in 
SALS one has to know the precise relation between pixels and scattering angles 
and this is troublesome when the distance z is changed to select a new particle 
diameter instrumental range. Also, and more important, SALS is plagued by 
stray light. To mitigate its effects, one has to rely on blank measurements to be 
subtracted from raw scattering data. The trouble is that stray light is worst at 
smaller angles, where the sensing elements are necessarily in small number and 
crowded close to the optical axis. With the present technique, on the contrary, 
all the pixels are used in calculating the correlation function for any value of 
the displacement r and this allows more accurate stray light subtraction; the 
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Figure 2.4: Optical layout for ENFS. The main transmitted beam is not blocked 
by the stop in the focal plane. Both the scattered and the transmitted beams 
are sent to the CCD. 

algorithms to subtract the stray light will be described in Chapter [31 

The results of the measurements on some colloid samples are presented in 
Chapter [71 The ONFS technique in the present form has only one tight require- 
ment, namely the clean disposal of the transmitted beam that requires accurate 
focusing and a proper diffraction limited beam stop. It is both conceptually 
and in practice very simple, and it capitalizes on the high statistical accuracy 
permitted by the large number of pixels of a CCD and by the good handling 
capabilities of PCs. 

It became soon appcarcnt that the main problem with ONFS comes from 
the poor statistical quality of the calculated / (q) . In Chapter [3 we will show 
that the statistical quality increases only as the fourth root of the number of 
processed images. We experimented a different optical setup (ENFS), drawn in 
Fig. 12.41 In ENFS, there's no beam stop: the main beam is let interfere with 
the scattered light. This is basically an heterodyne version of NFS, thus we call 
it hEterodyne Near Field Speckles (ENFS). 

Basically, ONFS data processing consists in evaluating the field correlation 
function Ce (r) by using Siegcrt relation (|2.2[) , then evaluating / (q) by applying 
the inverse Fourier transform to (|2.3p . In ENFS, we measure the interference 
between the speckle field of ONFS with the much more intense transmitted 
beam. We directly measure a quantity linearly related to the field. The intensity 
correlation function of an ENFS image equals Ce (r), provided that all the 
conditions needed by ONFS are met, that is, if the field is circular gaussian. We 
thus obtain C e (r) without the data inversion needed to apply Siegert relation, 
and this greatly enhances the statistical accuracy of the results. 

In Chapter [71 we show a comparison between data taken with ONFS and 
ENFS; data taken with ENFS are evidently much less noisy. The quality is 
comparable with the SALS one. This good quality allowed to try a Mic-based 
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Figure 2.5: Optical layout for SNFS. Part of the transmitted beam is blocked by 
a blade in the focal plane, along with half of the scattered light. Only one half 
of the scattered light is sent to the CCD, along with a part of the transmitted 
beam. 

inversion algorithm, to obtain an histogram of the distribution of the diameters 
of some colloidal samples; the measurements are shown in Chapter [5] 

Both ONFS and ENFS are quite sample wasting techniques. They require 
a sample much bigger than the statistical quality needs. For example, con- 
sider a non-equilibrium fluctuation measurement in a free diffusion experiment 
[13]. The biggest fluctuations we want to measure are about 0.5mm. A good 
statistical sample should be so big to contain some hundred of the biggest fluc- 
tuations: it can be a square with a 5mm side. This is enough for SALS, but not 
for ONFS nor ENFS. In Chapt. [3] we will show that, if we want to cover two 
decades in wavevectors, we must use a sample with side 5mm • 10^. To cover 
two decades, we need a half a meter wide cell, while with SALS we can work 
with a half a centimeter wide cell! This is not a difficulty for particle sizing 
applications, but can become a serious problem when we want to analyze many 
lenghtscales, since NFS is particularly suited for big objects. This problem is 
essentially due to the fact that big objects need long values of z in order that 
their scattered field is gaussian; on the other hand, we need a big sample, so 
that the sensor collect the light scattered at high angles by small particles. This 
fact is quite unusual, since in general big objects are good subjects for classical 
microscopy techniques. The difficulty can be easily circumvented, introducing 
a new instrumental setup, called SNFS: see Fig. 12.51 This setup closely mirrors 
the Schlieren setup: a blade stops half focused trasmitted beam, along with one 
half of the scattered light. In the case of a two dimensional sample, we can 
create an image of it on the CCD sensor, as in classical Schlieren imaging tech- 
nique. But in the case of three dimensional samples, the smaller objects create 
a speckle field, while the bigger ones are completely resolved: this technique 
passes continuously from the deterministic image formation of the big objects, 
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like in microscopy techniques, to the analysis of stochastical interference pat- 
terns of NFS, without the need to know at what lengthscale the passage takes 
place. We will show that SNFS doesn't require any condition on the position 
of the focal plane: it can be made as near to the sample as we need; moreover, 
the dimensions of the sample can be quite small, as in SALS. 

SNFS requires an additional element with respect to ENFS, the blade, but it 
allows easy measurements on many lengthscales, on big objects. We used such 
a technique to measure the power spectrum of non-equilibrium fluctuations in 
a free diffusion experiment, described in Chapter [HI thus showing that this 
technique can be applied to researches in fundamental physics. 
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Chapter 3 

Theory. 



In this chapter we describe quantitatively the relations between the scattered 
intensity and the properties of the field on a given plane. We will derive the 
transfer function of various techniques based on the light scattering. We will 
show that, under suitable conditions, the field is gaussian. In this case, we 
will derive the relation between the power spectra of the NFS images and the 
scattered intensity. 

3.1 Scattered intensity and field power spectrum. 

From Maxwell equations, we can derive the wave equation for a transverse 
component of the electric field in the vacuum |14j : 



2 



d 

-^E{x,z,t)=c^ 



,2 



d 

—^E{x,z,t) + ViE {x,z,t) 



(3.1) 



where V| is the Laplacian oeprator, with respect to the horizontal coordinate 
X. Since we are working with a laser, and we consider only elastic scattering, 
the only temporal frequency involved is kc: 

E{x,z,t) = E{x,z)e^^'^^^ (3.2) 



where k is the wave vector. Eq. (|3.ip becomes: 

— E (.f , z) + ViE (f , z) + eE {x, z) = (3.3) 
oz'^ 

Wc define Ez (q) as the Fourier transform of E {x, z) with respect to x: 

^Eziq) = -{k^-q')EAq) (3.4) 



The solution is: 



Ez{q)=E,{q)e'V}^^^ (3.5) 
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In order that this sohition exists, a condition must be fulfilled: 

< k^, (3.6) 

This condition is alwais met if we consider only propagating waves. 

The quantity Ez (q) is closely related to the intensity of the light crossing 
the plane z = cost. Each two-dimensional Fourier mode of amplitude Ez (cf), on 
a given z, and wavevector q is generated by a three-dimensional plane wave of 
wavevector [q^, qy, kz], where the only value of kz is obtained by imposing that 
the wavevector of any plane wave has length k: 

kz (q) - (3.7) 

Given the values of the two-dimensional Fourier modes Ez (q) on a given z ~ 
0, we can evaluate E(x,z) for each x and z, by using Eq. (|3.5|) and (|3.7p . 
Expressing it by its three-dimensional Fourier transform: 

E (q, kz) = 27r£;o (q) S [kz - fc, {q)] (3.8) 

Each three-dimensional component with amplitude E (g, kz) of the electric field 
represents a plane wave travelling in a different direction. We define Se {q, qz), 
the two-dimensional power spectrum of E (x, z): 

{Ez {q) E: (q')) = S {q-- q') Se (q) (3.9) 

Light intensity, for each scattering direction, can be defined on the basis of 
E{q,kz): 

{E {q, kz)E* {(f, k'J) = AttH {q~ cf) 6 {kz ~ k'J 6 [kz ~ kz {q)] I {q, qz) (3.10) 

where / {q, qz ) has been expressed in terms of the transferred wavevector [q^ ,qy,qz] = 
[kx,ky,kz]- [0,0, k]. FromEq. ([SJ]) : 

Qz {q) = Vk^ -q^-k (3.11) 

Substituting E {q, kz) of Eq. (|3.8p in Eq. (|3.10|) . and comparing the result 
with Eq. (|3.9[) , we can relate the scattered intensity / {q, qz ) to the power 
spectrum of the field Se (q)- 

I{q,qz)^SE{q) (3.12) 

From Eq. (j3.12p we notice that / {q) can be measured by evaluating Ez {q) 
on any z. 

If the sample is isotropic, / {q, qz) depends only on Q = \qx, qy, qz\'- 



Q(g) = x/2fcyi-^l-(|) (3.13) 
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Figure 3.1: Relation between (/ and Q. Geometrical interpretation of Eq. (|3.13p 
In this case, Eq. (|3.12p can be written in terms on q and Q: 



The geometrical meaning of Eq. (|3.13p is explained in Fig. 13.11 For q ^ k, Eq. 
(|3.13p can be approximated by Q (q) = q. Moreover, if Raylcigh Gans approxi- 
mation holds, / (q) represents the power spectrum of the refraction index of the 
sample. From these two considerations, we obtain the result that, for scattering 
on small angles and under Rayleigh Gans condition, the two dimensional corre- 
lation function of the electric field is proportional to the correlation function of 
the light path through the sample. 

When performing a far field, small angle scattering measurement, the scat- 
tered beams are focused on a screen. In suitable units, each point of the screen 
has a coordinate q. For small values of the wave vector, q approximates Q', the 
transferred wavevector. The exact relation is: 



Equation (|3.13p can be used to correct the results of a Near Field Speckles 
measurement. Figures 13.11 and 13.21 show the geometrical meaning of equations 
(|3.13[) and (|3.15[) . For small values of q, that is q/k ^ 1, the two equations 
can be approximated with Q = Q' ~ q', the error due to this approximation is 
shown in Fig. 13.31 it's quite small, and it can often be neglected. 



/ [Q {q)] = Se (q) 



(3.14) 



Q' (q) = V2k 1 



1 



(3.15) 
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Figure 3.2: Relation between the coordinate g on a screen, in a far field exper- 
iment, and the transferred wave vector Q' . Geometrical interpretation of Eq. 




0.01k 0.10k 1.00k 

q 

Figure 3.3: Relative error obtained neglecting the non linearity of the ralation 
between the sample wave vector and the near field wave vector. The graph is 
obtained from Eqs. (1X1^ and 
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3.2 Scattering from a thin sample. 



Let us consider a thin sample, with a non homogeneous refraetion index, and a 
hght plane wave, moving in the direction of the z axis. For z = 0, at the surface 
of the sample, the field will be: 

E{x,z = ^) = E^^^^^^)^ (3.16) 

where 81 (x) is the difference between the light path, the integral of the refraction 
index along z, for a given point x, and its mean value over the whole sample. 
If 51 is small compared to the light wavelength, we can consider a first order 
developement: 

E{x,z ^0) Eo[l + i6l{x)k] (3.17) 

Neglecting the higher order terms means that we are neglecting higher order 
diffracted beams than the first. 

Using Eq. (j3.5p we can find the field for every value of z: 

E,ix) = Eoiz) + 6E,{x) (3.18) 

where 

SE, (q) = ikEo (z) SI [q) e (^^^ ^^^^^^ 
is the scattered field, and 

E^ [z) = E^e'^^ (3.20) 



3.3 Image forming techniques 

Microscopy, in its basic form, consists in forming an image of a plane on a 
device which measures light intensity, such as a photographic film or a CCD 
sensor. Generally it is used to obtain informations about the intensity of the 
transmitted light, for example, in the case of an organic tissue, treated by some 
dye. 

If a microscope objective forms the image of a plane on the CCD, the image 
is given by the interference of the transmitted and the scattered beams. For such 
an image, the signal can be defined as the difference of the measured intensity 
I {x) and the transmitted beam intensity /q, divided by Iq. We will call this 
signal ishadowgraph, for rcasons that will be clear later. We consider the case in 
which the scattered beams arc much less intense than the transmitted one. At 
the first order in SE, the signal is: 

ishadowgraph (x) = — = [Eq5E* (x)] (3.21) 

JO Jo 

and the Fourier transform is: 

ishadowgraph {q) = ^ [EqSE* {-q) + E*5E {q)] (3.22) 
Jo 
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If the sample is transparent, and we send a plane wave through it, Eq. (|3.19p 
tells us the value of the scattered field: 

ishadowgraph z) = 2kSl (cf) sui (^k - \/k^~(f^ z (3.23) 

In order to obtain the previous result, bl (x) has been considered real, so that 

For z = 0, that is, if the thin sample is in the focal plane, ishadowgraph {q, z) — 
0. The intensity is completely uniform, and bears no informations on the sample. 

Many techniques has been developed in order to make phase modulations 
evident: among them, holography and interferometry. A well known way to 
make phase modulations evident is the phase contrast microscopy. Basically, 
this technique consists in changing the phase of the transmitted beam by 7r/2. 
At the first order in SE: 

iphase contrast (x) = — = ^3 [EqSE (f)] (3-24) 

Jo Jo 

and the Fourier transform is: 

iphase cor^trast io) = ^ [-^^E {-q) + E^SE* {q)] (3.25) 

Jo 

Using Eq. (|3.19p to evaluate the scattered field: 

iphasecontrastiq,z) =2k6l{q) cos (^k - y/k^~~(f^ Z (3.26) 

For z = 0, that is, with the sample in the focal plane: 

iphase contrast {q, Z = 0) = 2kSl (q) (3.27) 

Another way to make phase modulations evident is the so called dark field 
technique. It consists in stopping the transmitted beam. This is accomplished 
by focusing the transmitted and the scattered beams by a lens, and by removing 
the transmitted beam by some kind of refiecting or absorbing object. This is 
an homodyne technique; the signal must be defined as the ratio between the 
measured intensity / (x) and the intensity of the transmited beam Iq. Since we 
have only 5E: 

I{x) \5E{x)\' 

^dark field \X) = — — — {6.Zii) 

Jo Jo 

Equation (|3.19p . for z = 0, gives: 

idark field {x, z = 0) = k^ (x) (3.29) 

In Fourier space: 

idark f^elf (<f, Z = 0) - T^^" / i^) ^1 (<?- dcf (3.30) 

{2tt) J 
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The Schlieren technique consists in focusing the beams from the sample by 
a lens; in the focal plane, a blade stops half of the transmitted beam, along with 
the beams scattered in one half plane. At the first order in SE, the signal is 
again, like in shadowgraph: 



(x) 



I{x)-Io 



In 

lo 



(3.31) 



where the field SE is the scattered field, without one half plane in Fourier space: 



6Eiq} = 



5E{q) q-n<0 
SE{q) q-n>0 ' 



(3.32) 



n is the vector ortogonal to the direction of the blade, Iq = Iq/2 and Eq = 
Eq/^/2 are the intensity and the field of the trasmitted beam, after the blade. 
The Fourier transform is: 



ISchli 



Using Eq. (^3^ : 



'^Schli 
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EoSE* i^q) + E*SE (q) 



^Eo6E*{-cf) q-n<Q 
^E*5E{q) q-n>0 



(3.33) 



(3.34) 



Using Eq. ^JQ: 



^Schli 



V2ikSl (q) e 
V2ikSl (q) e 



— q'^ — k] z 



\/k'^ — q^ — k \ z 



q-n<0 
q-n>0 



(3.35) 



In order to obtain the previous result, 51 (x) has been considered real, so that 
61* (—q) = 6l{q). For z = 0, that is, if the thin sample is in the focal plane: 



iSchlieren iq, Z ^ 0) = V2ikdl [cf) 



(3.36) 



The factor i means that all the Fourier components of SI (q) undergo a rotation 
of 7r/2: a sine-like light path modulation gives a cosine-like intensity modulation. 



3.4 Misfocused microscopy and shadowgraph. 

Equations (|3.23p , (|3.26p and (|3.35p allow to evaluate the evolution of the signals 
as z, the misfocusing, is increased, for a simple microscope objective, for phase 
contrast and dark field. It should be noted that we are dealing with images 
formed by laser light: as z increases, it is possible to recover the original shape of 
the observed objects; this is completely different from a white light microscope, 
in which the misfocusing simply smears the images. 
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Equation (|3.30p has not been extended for z 7^ 0; in all the above mentioned 
techniques, however, the variation of z strongly influences the relation between 
the light path {q) and the signal i (q). This is generally a defect: thick objects, 
or even thin objects dispersed in a thick volume, are difflcult to be analysed. 
In general, all the above mentioned techniques are applied to samples that are 
thin, and in the focal plane. No improvement is obtained by misfocusing. 

A well known exception is shadograph. Shadowgraph technique consists in 
sending a plane wave onto a sample, and observing the intensity modulations 
generated by the sample on a plane placed at a distance z from the sample. 
Using Eq. p.23p . we can derive the transfer function T shadow graph (9, z) of the 
shadowgraph technique [TH [THl E] : 

Tshadowgraph z) = 2fc siu - B - q^^ z « 2fcsin (3-37) 

The approximation holds for q <^ k. The transfer function is defined as the 
ratio between the signal and the light path modulation amplitude: 

tiq,z)^T{q,z)Sl{q) (3.38) 

For z = 0, the transfer function vanishes; misfocusing is needed, and is a 
simple way to make phase modulations evident. 

Looking at Eq. p.26p we can notice that phase contrast transfer function, 
as a function of z, has a cosinusoidal behaviour: 

Tphase contrast {q, z) = 2fc COS (^k - Z W 2fc COS (3.39) 

When considering opaque objects, with no phase modulations, the transfer func- 
tions are exchanged: cosinusoidal for shadowgraph, sinusoidal for phase con- 
trast. 

The shadowgraph image is created by the interference between every scat- 
tered beam and the transmitted beam; it's always possible, in principle, to find 
the value of 51^ in every point, simply by a deconvolution. Shadowgraph al- 
lows the measurement of one component of the field, which, in turns, is the 
convolution of the light path with a particular function. The absolute intensity 
modulation of the shadowgraph image is proportional to the mean intensity and 
to the light path modulation; the constant of proportionality is the the transfer 
function. The transfer function vanishes for some wave vectors, but has maxima 
for other ones. At the maxima, the sensibility equals the sensibility of phase 
contrast and Schlieren techniques. 

Now we evaluate the effect of misfocusing on a dark field microscope. We 
obtain an image of a plane a distance z from the cell. Using Eq. I|3.19p . we 
can derive the relation between the light path and the measured intensity, at a 
given z: 



-I 

idark field (?, z) = " o fc^ / 51 ((f , z) 51 (q - (f, z) 6 

(27r) 



^Jk^-iq-q'f 



z 



(3.40) 
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This expression reduces to Eq. (|3.30p for z = 0. At this point, there's no 
appearent reason to use a misfocused dark field instead of a focused one. 

The knowledge of idark field (si,z), for every z, can give some informations 
about the spreading of the scattered light. For the scattering of a single particle, 
one can measure the intensity on planes with increasing values of z, and calculate 
how fast the light is diverging. This provides informations both on the position 
of the particle and on the scattered intensity. For a sample composed by a great 
number of particles, this cannot be done, and a different, statistical approach 
must be applied. 

For Schlieren technique, from Eq. (|3.35p : 



rp \ J V^ike V / q-n<0 , , 

J-SchHeren[q,Z) - / \ (3.41) 

V2ike V ^ ) ^ q-n>0 

By evaluating the square modulus of the transfer function, we obtain: 

\TschHeren{q,z)f =2k^ (3.42) 

This means that the power spectrum of the electric field is proportional to the 
power spectrum of the light path, without any dependence on spatial wavelength 
and misfocusing. 



3.5 Scattering measurements by microscopy tech- 
niques 

By using Eq. (j3.12p . we can, in principle, evaluate the scattered intensities by 
measuring the field on a given plane. The above described microscopy techniques 
allow the determination of some functions of the electric field. Now, we assume 
that the scattered electric field is generated by a thin sample. From Eq. (|3.19p . 
we derive a property of the electric field: 



6E i-q) = -SE* (q) e 



2^{V^-k)z 



(3.43) 



For shadowgraph, we use Eq. (|3.22p . in order to evaluate the power spectrum 
Si of i (x) : 

4 



Jo 

Using Eq. p.l2p . we obtain: 



k — ^Jk"^ — q^ ] z 



Se (q) 



(3.44) 



S,{q) = sm^ [k-^k^-q 



i[Q (q)] 



(3.45) 



In order to use shadowgraph to evaluate the scattered intensity, the sample 
must be out of the focal plane. This technique has some disadvantages: some 
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wavevectors cannot be seen, since the transfer function vanishes; moreover, if 
the sample is thick we cannot define a z: the oscillations of the transfer function 
are smeared, but it could be hard to know quantitatively how much. 

For phase contrast, wc use Eq. (|3.25|) . in order to evaluate the power spec- 
trum Si of i (x): 



Jo 

Using Eq. (|3.12p . we obtain: 



(fc - ^k^ - z] Se (q) 



(3.46) 



Si {q) = -r 
Jo 



k-^/k^ 



I z 



liQiq)) 



(3.47) 



The disadvantages of shadowgraph, as a scattering measurement technique, are 
also found in phase contrast. Phase contrast has a flat transfer function only if 
z = 0: it's an ideal technique for thin samples. 

Dark field doesn't allow to recover any information about the phase of the 
field. It is not suited to make scattered intensity measurements; a statistical 
approach, described in the following sections, will give interesting results. 

We will describe in Section [3.1Il the application of Schlieren technique to the 
measurement of the scattered intensity. 



3.6 Gaussian field generated by the sum of many 
patterns. 

Consider a monodispcrsc colloid, and the near field scattered light. The field 
can be decomposed into the sum of the waves coming from the different elements 
of the colloid. Thus at every distance from the colloid, the field will be given by 
the sum of many patterns, each randomly placed. An example of this is given 
in Fig. 13.41 and 13.51 Figure 13.41 shows the pattern: the intensity is linearily 
dependent on the field / (x), which we can consider given by the wave scattered 
by a particle of the colloid. Figure 13.51 shows the sum of the patterns, the 
function p{x): 

N 

p{x) = Y,f{x-x,) (3.48) 

1=1 

where the N randomly distributed in a surface of measure S. 

Now wc will evaluate the P-point correlation functions of the sum of many 
patterns, p{x). The results will be obteined, first, for fixed particle density 
J\f = N/S, and N 00; then we will show that, in a suitable limit, the 
P-point correlation functions, for P > 4, can be expressed in terms of two- 
point correlation function, corresponding to the Wick formula. In other words, 
every connected part of the correlation function developement vanishes: the 
field becomes gaussian. As a matter of fact, we will prove an extension of the 
well known central limit theorem. 
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Figure 3.4: Example of the field generated by a particle of the colloid. 
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In the following, we will consider only functions with a vanishing average 
value, the other cases being easily obtained from this one. This simplifies the 
problem, since every odd-P-point correlation function will vanish. 

The P-point correlation function of p {x) is: 

C(Ax2, . . . , Afp) = (p(.f)p(x + Ax2)...p(x + Axp)) = ^ 



gN+l 

N 



V //(f-£.)/(^_x, + Af.).../(f-^. + Afp)dfdf,...df. 
Js 

(3.49) 



i,j,k,--- = l 



The value of the integral docs not depend on all the values of the indices 
k, . . . , but only on which of them are equal; the sum involves NP terms, 
but many of them are equal. For example, for P = 4, the term with i = l,j = 
2,/c = 3, / = 4 is equal to the one with i = 2,j ~ 5,k ~ 7,1 ~ 9, but it is 
different from the one with i = 1, j = 1, fc = 3, ^ = 4. The problem is thus to 
determine in how many ways we can obtain a given configuration. 

The calculation can be made more easy using graphs. For evaluating a P- 
point correlation function, we draw P points on a graph, each one corresponding 
to one of the points of the correlation function, 0, A.T2j • ■ • , Axp. Then, we 
group the points, so that every set contains an even number of points 0. Each 
configuration corresponds to many values of the indices i,j,k,...; the number 
of them is the multiplicity of the graph. Every set corresponds to an operation 
of integration on a different Xi. We call G the number of sets; we will have only 
G integration variables, being the integrand independent on the other N+l — G 
variables. The integration on these variables gives a factor S^'^^~'~^. Moreover, 
every integration corresponds to the evaluation of the correlation function C of 
the single pattern / (x) : 

C iAx2, Axp) = ijf{x)f {x + Ax2) ...f (x + Axp) dx (3.50) 

The multiplicity of the graph depends on G; its value is iV (iV - 1) (iV - 2) . . . (iV - G + 1) 
For N —> oo, we can consider only the leading term N'~^. 

We can thus describe the rules for evaluating the correlation functions, as 
the sum of all the graphs. The value of every graph is the product of the factors 
given by each set. The factor is the product of J\f and the correlation function 
G, which correlates all the points in the set. 

For example, we evaluate the two-point correlation function (Tab. 13. ip and 
the four-point correlation function (Tab. 13. 2[) . 

For A/" — > cxD, the leading term in the developement of the correlation func- 
tion is the one with the higher power of Af: it is the one with the higher number 
of sets. Since sets with an odd number of elements have a vanishing contri- 
bution, the greatest number of sets can be obtained only by making sets of 



^ Sets with odd number of points will give vanishing contributions if / (x) has a vanishing 
mean value. 
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10 Afil C(Axi) AAC(Afi) 

Table 3.1: Evaluation of the two-point correlation function. 





Afg 





Axi 



I Ax2 Ajfjl 



|0 Afil 




C(Afi,Ax2,Af3) . 
AfC (Afi, Af2, Af3) + 
AA2C(Afl)C(Ax2- Af3)- 
A/'2(7(Af2)C•(Afl - Axa)- 



Af^C (Afg) C (Af 1 - Afs) 



^a;2X^2;3\ 

; XAf 1 

Table 3.2: Evaluation of the four-point correlation function. 
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two points. This means that only two point correlation functions of the single 
pattern C (Aaf) contribute to any correlation function C (Aafi, . . . ) of the sum. 

Since M is dimensional, it is not possible to state if it is small or great. This 
means that wc cannot define, in general, a value of M so great that the field 
becomes gaussian. The following heuristic considerations will show that the 
field is gaussian if MA ^ 1, where A is the area of one pattern, at least if we 
can define it in some ways. Consider a pattern / {x) ~ axA (x). The P-point 
correlation function of the pattern /(x), evaluated in Ax = 0, has the value 
(7 (0, ... ) a^A. Every graph will have a factor and a factor A^, where G is 
the number of sets in the graph. So the factor J\f alwais appears multiplied by A. 
By imposing MA ^ 1, we obtain that the only contributions to the correlation 
function of the sum of patterns comes from the two point correlation function 
of the single pattern: all the Wick formulas are valid. In order that JVA ^ 1, 
the mean number of scattering particles inside each area A must be large: many 
pattern must overlap, in each point. 



3.7 Siegert relation for the near field speckles. 

Each scatterer of the sample generates a diffraction pattern which, at least in its 
far field, becomes larger and larger linearly, as the distance z from the screen and 
the sample is made longer. So, for z longer than a given distance, many diffrac- 
tion patterns overlap: the Wick formulas should become valid. Unfortunately, 
the considerations of the previous section cannot be applied directly, since the 
area the diffraction pattern has not been already defined in a quantitative way. 

Now we will prove the Wick formula for the case of Siegert relation, by using 
the formalism developed in the previous section. 

From the results of the previous section, we can evaluate the intensity cor- 
relation function of the sum of the patterns, at a given z: 



\EAx)f\EAx + Ax)f 



(x) E^ (f + Ax) 



dx 



j (y + A 



\x) 



dxdy ' 



Af^ / E^ {x) El {x + Ax) El (y) E, {y + Ax) dxdy ^ 



/ E, (f) E, {x + Ax) El {y) El [y + Ax) dxdy 



(3.51) 



where E^ (x) is the field from a single scatterer, at a distance z. Assuming that 
the phases are random, and considering that the correlation function of any field 
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does not change with z: 



{h (x) h {x + Ax)) = N h {x) h {x + Ax) dx + 



(3.52) 



Jo {x) dx 



Eo (x) Eq {x + Ax) dx 



The first term on the right hand side depends on z: it is the intensity correlation 
function of the diffraction pattern. Since the diffraction pattern becomes larger 
as z increases, while the total intensity keeps its value, the intensity correlation 
function of the diffraction pattern decreases, and vanishes as 2; ^ 00. 

In order that Siegert relation holds, for a finite value of z, we must impose 
that the term with the four point correlation function is negligible compared to 
the two point ones: 



A/" / (x) (x + Ax) dx < 7V^ 



7o (x) dx 



(3.53) 



The first term can be substituted by its higher value, the one with Ax = 0: 

n 2 



Iz {x) dx 



> 1 



If (a-)df 



The fraction represents the area A covered by the diffraction pattern: 

AAyl> 1 



(3.54) 



(3.55) 



In order that Siegert relation holds, we need that many particles scatter light 
inside a single diffraction pattern, that is, any point of the screen must be hit 
by light coming from many particles. This can be obtained without changing 
Af, but simply increasing z, thus increasing A. 

It should be noted that the validity of Vick formulas for a given z does not 
mean that the field is completely gaussian. For example, we have shown that 
Cj (Ax) for 2 — + 00, where Cj (Ax) is the intensity correlation function 
of the diffraction pattern, defined as J E [x] 1^ {x + Ax) dx. But this does not 
imply any uniform convergence. Its integral, / C| (Ax) dAx, for example, is 
a constant, and does not vanishes as z 00. This means that we can build 
suitable linear operators, acting on the field, yelding quantities which do not 
have a gaussian distribution. A dramatic example can be obtained considering 
the scattering from a two dimensional screen, with many holes of a given shape. 
As z — > 00, the field meets the Vick formulas ever better. But it is alwais 
possible to analyze an area, bigger than the diffraction pattern of each hole, 
and to recover the shape of the holes. This can be done by deconvolving the 
field by a suitable function: it's the operation made by a lens, which creates an 
image of the holes. The deconvolution gives any information about the sample. 
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including the fourth order correlations: the deconvolved field is not gaussian. 
The gaussianity is only local: once we defined an area, corresponding to the 
aperture of a lens, there's a distance beyond which we are not able to recover 
the shape of each hole, and so informations on higher order correlation functions 
than second order ones are lost. 

We can conclude that Eq. p.55p implies only a local gaussianity; gaussianity 
is valid only when considering points inside an area small compered with the 
diffraction pattern of each scatterer. On the other hand, the knowledge of 
the field on a whole plane allows to recover any information on the correlation 
function of any order. 

3.8 Vanishing of the {EE) correlations. 

In Eq. (|3.5ip we neglected the terms like / E [x) E {x + Ax) dx. Such terms 
should give contributions like {E (a?) E (x + Ax)) in the Vick formulas; wc ne- 
glected them assuming that the phases are random. In this section we will 
analyze the conditions under which this happens. 

In general, (EE) correlations are not negligible. For example, we can con- 
sider an opaque screen, with a transmission coefficient dependent on the posi- 
tion, with gaussian distribution. We know that Vick formulas hold. Now, we 
send a beam through it, and measure the outgoing intensity correlation function, 
immediately after the screen. Since all the points are in phase, we can assume 
that the field is real. The Vick formula for the intensity correlation function 
states that C/ {Ax) = (/)^ -I- 2 \Ce (Aa;)|^, due to the not neghgible contribu- 
tion of the term (EE), which becomes equal to (EE*). Another example can 
be found in the theory of shadowgraph: the term (EE) is responsible for the 
oscillations of the transfer function defined in Eq. ()3.37|1 . 

Now we derive the equations giving the evolution of (EE) as z increases. 
We define: 

F, {Ax) = J (^) {x + Ax) dx (3.56) 

The Fourier transform of F {Ax) is: 

{q) = E, {q) E, {-q) (3.57) 

We can notice that F {q) is the power spectrum if E {^q) is the complex conju- 
gate oi E {q), that is if E {x) is real. By using Eq. (|3.5p . we obtain the evolution 
ofF,{q): 

F. {q) = e'^'V^^^Eo {q) E, {^q) (3.58) 
This gives the evolution equation of {q): 

{q) = e2Vfc^ - (g) (3.59) 

Comparing this equation with Eq. (|3.5p . we see that F^ {q) evolves like the 
electric field Ez {q), but two times faster than it, as z increases. 
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The root mean square amplitude of (q) is a conserved quantity; since 
Fz (x) gets larger and larger as z increases, its amplitude must decrease like 1/z. 
We can thus define a condition which is enough to ensure that the terms (EE) 
are negligible: the diffraction pattern must be much larger than the correlation 
lenght. This is implied by Eq. (|3.55|) . The gaussianity condition expressed 
by Eq. (j3.55p is met if many diffraction patterns overlap in every point. This 
implies that the diffraction pattern of each object must be much larger than the 
object itself, and than its correlation function, at least if the objects themselves 
do not overlap. 

Some difficulties arise when we consider the power spectrum, or the Fourier 
transform of (EE) terms. As we already explained, the root mean square value 
of (EE) does not depend on z. A Fourier transform, made over a whole plane 
at a given z, could be divergent, for some values of q, as z increases. For 
example, we consider a Fourier transform made on a given area S, and we 
evaluate its mode with wavelength 0, that is, the integral of Fz (x) over S. 
It is proportional to S/z. We can consider a square area S, of side 2'K/qi, 
where qi is the wavevector of the longest wavelenght Fourier mode of the square 
S. So Fz {q~0) cx l/qfz. Once we selected a qi, the lowest wavevector we 
will consider, in order that (EE) is negligible with respect to a given value, 
independent on z, we must impose a z cx 

A more quantitative result can be obtained by considering the evaluation of 
the Fourier transform on S as the evaluation of the Fourier transform on the 
whole plane, followed by the convolution with the Fourier transform of xs i^)- 
This is equivalent to considering the discretization of the allowed wavelengths, 
due to a finite area S. Near a given value of q, the exponential term in Eq. (|3.59p 
makes an oscillation in about k/qz. The discretized intervals are spaced by qi: 
if k/qz <C qi the oscillations are avereged and vanish. In general, the oscillations 
will be more visible for small values of q. In order that the oscillations are never 
visible, k/qiz <^ qi. once we have selected gi, that is the side of S, we must 
provide that: 

z>4- (3-60) 

In shadowgraph technique, Eq. (|3.60p means that the oscillations of the 
transfer function are so fast that they cannot be resolved by the sensor, and are 
thus averaged. 

Equation (|3.60p has a geometrical interpretation. The vanishing of {EE) 
can be expressed in terms of Fourier modes: 

{E{q)E{-q))={) (3.61) 

The beams, scattered by a modulation with wavevector g, hit the sensor at 
an angle q/k. Every modulation with wavevector q scatters at two simmetric 
angles; the resulting modulation on the sensor is thus given by light coming 
from two different regions, of area 5, whose distance is about Iqzjk. If the 
distance is longer than 2iT/qi, the side of S, the regions do not overlap: see 
Fig. 13.61 Equation (j3.60p states that we must provide that the regions do not 
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Figure 3.6: Description of the condition of non overlapping of the scattering 
regions. 

overlap. This means that light collected at symmetric angles is not correlated, 
as required by Eq. p.6ip . This condition ensures that the field is gaussian, only 
if the density of scatterers M is so that MS ^ 1 . 

For the scattering from a thin sample, the intensities of the beams scattered 
at two symmetric angles are equal, and the phases are defined. Since the angles 
are symmetric, the interference of the scattered beams with the much intense 
transmitted beam gives two interference patterns, sinusoidal modulations, with 
the same wavevector, and a given phase. Changing z, the two diffraction pat- 
terns change their phase; at some z they sums, and at other values they cancel 
out. This is the origin of the oscillations in the transfer function of shadow- 
graph technique, defined in Eq. (|3.37p . If the condition of Eq. (|3.60p is met, 
the phases of the beams scattered at symmetric directions is random: on aver- 
age, the transfer function is constant. 

The vanishing of the (EE) terms can be obtained also by increasing the 
thickness of the sample. When we pass from the two dimensional, Raman 
Nath scattering to the three dimensional, Bragg scattering, the correlations 
between the two beams scattered at the symmetric angles by a given sinusoidal 
modulation are not preserved. In shadowgraph language, the transfer function 
oscillations are washed out by superposing many layers, at different z. The 
thickness of the sample 5z must meet the condition Sz > k/ q^. 

3.9 Homodyne near field speckles. 

This technique has been presented very recently [U [2]. The device for the 
measurement of the Near Field Speckles is, basically, a misfocused dark field 
microscope. The transmitted beam is removed, and the image is due only to 
the light scattered from the sample. 

Some of the parameters of the system must be selected: the distance z from 
the sample to the focal plane of the objective, or from the CCD, if the objective 
is missing; the diameter D of the sample and of the incident beam; the superficial 
particle density of the sample. The parameters must be selected on the basis of 
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the required wavevector range [gmm, 9maa:]- The ratio Qmax/qmin cannot exceed 
two decades, due to the finite size and discretization introduced by the CCD 
sensor. The range [qmim Qmax] is generally selected in order to cover interesting 
wavelenghts of the sample: for example, from one tenth to ten diameters of 
the particles, in the case of the scattering from a monodisperse colloid. Three 
conditions must be fulfilled. 



1. The misfocusing z must be selected in order to meet the condition ex- 
pressed by Eq. (|3.60|1 . with qi = qmin, in order that the correlations {EE) 
vanish: 

z > (3.62) 

This condition is stronger than z 3> k/ (qtypqmin), where 2Tr/qtyp is the 
typical diameter of the particles. Since iir/qmin is the side L of the images 
we take, Eq. p.62p implies that qtypz/k 3> L: the diffraction pattern 
of each particle covers a surface much bigger than the observed one, as 
required in order that the Fourier transform of the field can be considered 
gaussian. 

2. In order that the field is gaussian, many diffraction patterns must overlap, 
in each point. Under the condition of Eq. p.62p . in order to fulfill Eq. 
(|3.55p . we must only provide that there are many particles in the surface 
S covered by an image: 

MS > 1 (3.63) 

. For particles suspended in a three dimensional volume, we can define the 
superficial particle density by multipling the volumetric particle density 
by the thickness of the cell. In order to fulfill condition (|3.63p . we can 
increase the volumetric particle density, or increase the thickness of the 
sample. Care must be teken in order to avoid multiple scattering. 

3. The images we take must collect light scattered at any angle by the sample. 
The highest wavevector we want to measure is qmax] in order that the 
sensor collect light scattered by that wavevector, coming from any area of 
the sample, its diameter D must satisfy: 

D » (3.64) 

This condition ensures that the sensor cannot see the sample boundaries: 
the sample can be considered as infinite. If the diameter is much less than 
the one imposed by Eq. p.64p . the speckles are governed by the classical. 
Van Cittert and Zernike theorem. 



Under Eq. dSH]), (piBg]) and (PTM)) . Siegert relation holds: 

Ci{/\x) = (/(a-)/(x + Af)) = {if + \Ce{/^x)\^ (3.65) 
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where the (•) is the mean over x and Ce (Aa?) is the field correlation function. 
The intensity we measure in a point is not directly connected with any physi- 
cal part of the sample: each speckle is generate by the superposition of many 
diffraction patterns. 

The measurement of the intensity allows to recover the modulus of the field 
correlation function through Eq. (|3.65p . Since Ce (Ax) is the Fourier transform 
of the power spectrum, which is symmetric and real, Ce (Ax) is also real. More- 
over, if we think Ce (Ax) is alwais positive, we can calculate it by extracting a 
square root. Then, we Fourier transform Ce (Ax), thus obtaining Se (q)- 



Se (q) = ^ 



^Jci{x)-{iy 



(q) 



Using Eq. (|3.19p . we obtain, for a thin sample: 



^Ci (x) - {If 



(q) 



Using Eq. (|3.14p . we obtain the scattered intensity 

I [Q (<?)] = ^ 



Ciix)-{iy 



(q) 



(3.66) 



(3.67) 



(3.68) 



The results do not depend on z. The misfocusing z must be enough, in order 
that the field is gaussian, but its value does not affect the results. 

The extraction of the square root of the difference between two experimental 
data is a dangerous operation, since the difference could be negative. In gen- 
eral, as any other inversion of experimental data, it involves an increase and a 
distorsion of noise. Chapter [7] provides a detailed description of this problem. 
This kind of problems are avoided by using ENFS or SNFS, described in Sects. 
[3T0l and [3TT] 

Wc must notice that Eqs. (p?62)) and ((3^ give D > qmax/qlmn- We 
can think qmax / q-mm as the wavevector dynamic range we want to measure; 
hopefully it can be about one hundred for spatial measurements. On the other 
hand, I-k jq^nm is of the order of some length of the particles, for example ten 
times. For example, if we consider 10/im colloids, each image must cover about 
100/Ltm. In order to cover two decades in wavelength, we need a D about one 
hundred times wider: about 1cm. This is not a huge length; moreover, from the 
industrial point of view, there's no problem in making many acquisitions, with 
different magnifications and z, for every wavevector range, since each acquisition 
needs no accurate positioning. 

On the other hand, in some cases, for scientific purposes, D should be too 
wide. This is the case of measurements of non-equilibrium fluctuations, de- 
scribed in Chapt. O In this case, we want to evaluate power spectra on two 
decades in spatial frequencies. The dimension of the largest fluctuations is about 
one tenth of millimeter: a good statistical sample is about a millimeter large, 
and the whole sample must be two decades bigger: about 10cm. The building 
of a Soret cell, or a free diffusion cell, with such an big diameter can be avoided 
by using SNFS, which will be described in Section r3. Ill 
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3.10 Heterodyne near field speckles. 



We developed this technique very recently; the device has been patented [3l[4]. 
The setup for hEterodyne Near Field Speckles measurement is identical to ONFS 
one, but the beam stop is missing. The transmitted beam is not removed, and 
the image is due to the interference of the light scattered from the sample with 
the transmitted beam. 

Two parameters of the system must be selected: the distance z from the 
sample to the focal plane of the objective, or from the CCD, if the objective is 
missing and the diameter D of the sample and of the incident beam. In ENFS, 
the superficial particle density of the sample plays no role. The parameters must 
be selected on the basis of the required wavevector range [gmm, Qmax]- The ratio 
Imax/lmin cannot cxcccd two decades, due to the finite size and discretization 
introduced by the CCD sensor. The range [qmin, Qmax] is generally selected in 
order to cover interesting wavelenghts of the sample: for example, from one 
tenth to ten diameters of the particles, in the case of a monodisperse colloid. 
Two conditions must be fulfilled. 

1. The misfocusing z must be selected to meet Eq. p.62p . in order that the 
correlations (EE) vanish. This condition is stronger than z ^ k/ (qtypqmin), 
where 2Tr/qtyp is the typical diameter of the particles. Since 2TT/qmin is 
the side L of the images we take, Eq. p.62p imphes that qtypz/k » L: 
the diffraction pattern of each particle covers a surface much bigger than 
the observed one, as required in order that the Fourier transform of the 
field can be considered gaussian. 

2. The images we take must collect light scattered at any angle by the sample. 
The highest wavevector we want to measure is qmax] in order that the 
sensor collect light scattered by that wavevector, coming from any area 
of the sample, its diameter D must satisfy Eq. (|3.64p . This condition 
ensures that the sensor cannot see the sample boundaries: the cample 
can be considered as infinite. If the diameter is much less than the one 
imposed by Eq. (|3.64p . the speckles are governed by the classical. Van 
Cittert and Zernikc theorem. 

We can notice that the conditions expressed by Eq. (|3.62p and ()3.64|) must 
hold for both ONFS and ENFS. On the contrary, in ENFS no condition is 
imposed on the particle density, in analogy with (|3.63p . since the field does not 
need to be gaussian. 

In general, Eq. (|3.63p is fulfilled by the sample; in that case, the field is 
gaussian, and the ENFS image represents the interference of a gaussian field 
with a plane wave. The particle density can be so small that the field is not 
gaussian; this does not mean that the speckles we see represent real objects in 
the sample. Each speckle is due to the interference between light scattered by 
many different particles. 

Care must be taken in order to avoid multiple scattering. Since we want avoid 
multiple scattering, the scattered intensity is small compared to the transmitted 
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beam intensity: the second order effects in SE/Eq can be neglected, and Eq. 
[321] holds. 

C, (Af) = (Af) + [E*^ {SE {x) SE (f + A^))l (3.69) 

Under Eq. (|3.62p . the (EE) correlations vanish: 

a (Ax) ^ ^^CsE (Ax) (3.70) 
Jo 

This means that Si (q) = j^Se (?). The measurement of the intensity allows to 
recover the field power spectrum. 

Using Eq. (|3.19|) . we obtain, for a thin sample: 

2k^Ssi (q) = S, (q) (3.71) 

Under (|3.64p . Eq. (|3.14p holds. We obtain the scattered intensity: 

I[Q{q)]^^IoS^{q) (3.72) 

The results do not depend on z. The misfocusing z must be sufficent, in order 
that the correlations (EE) vanish, but its value does not affect the results. 

The considerations about the diameter D of the sample hold also for ENFS: 
Eq. (|3.62p and (|3.64p give D ^ Imax/qmin- This problem has been discussed 
in Section [3.9l The result is that, in some cases, the sample and the laser beam 
have to be extremely large. In that cases, SNFS can be used instead of ENFS: 
that technique will be described in Section r3. Ill 

Equation (j3.72p must be compared with Eq. (|3.45p . that holds for values of 
z much less than those imposed by Eq. (|3.62p . The oscillations in the sensibility 
of shadowgraph technique come from the non vanishing of {EE) correlations, 
essentially due to the phase relation of the beams scattered at symmetric angles 
by a thin sample. For example, the zeroes of the transfer function are due to 
the distructive interference of the symmetrically scattered beams. In ENFS, the 
phase relation is destroied, because the light that hit the sensor at symmetric 
angles comes from different regions. 



3.11 Schlieren-like near field speckles. 

For Schlieren technique, we can use Eq. (j3.34p . in order to evaluate the power 
spectrum Si of the signal i (x) : 

o / T;SE{-q) q-n<0 

I To {qJ q-'n>0 

We assume the sample is isotropic, so that / (q) depends only on |q|. Using Eq. 
(i3T4)l : 

I[Qm = ^hS^{q) (3.74) 
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Schlieren technique can alwais be applied to measure the scatterered intensity, 
no matter how long the misfocusing is. The sample can be thick or thin, in the 
focal plane or away from it: the result is never affected. 

Once the dimension of the image has been selected, in order to observe an 
interesting range of wavevectors, the diameter of the sample must fulfill the 
condition expressed by Eq. (|3.64[) . as in ONFS and ENFS, but in this case 
there's no limitation on z. The diameter will be, in general, sufficient to give a 
good statistical sample of the particles we are measuring; z will be as small as 
we can. 

In general, for a thick sample, some of the objects will be too small, or too far 
from the focal plane, to be completely resolved. But their presence will prodece 
a speckle field, analogous to that of NFS. We will call this technique Schliercn- 
like Near Fiels Speckles, since it behaves like a true Schlieren technique only for 
big objects in the focal plane, while for the other cases it allows to measure the 
statistical properties of a speckle field. 

Equation (|3.74p must be compared with Eq. (|3.45p . that holds for values 
of z much less than those imposed by Eq. (|3.62p . and without the blade. The 
oscillations in the sensibility of shadowgraph technique come from the non van- 
ishing of (EE) correlations, essentially due to the phase relation of the beams 
scattered at symmetric angles by a thin sample. In SNFS, the phase relation is 
destroied, because one of the beams scattered at symmetric angles is stopped. 

3.12 Why using NFS instead of classical scat- 
tering measurement? 

Classical, high angle light scattering works only for scattering angles higher than 
some degrees. Small angle light scattering (SALS) can measure light scattered 
by wavelengths of some microns; its main problem is the stray light. The new 
techniques allow an accurate subtraction of the stray light; in particular, for 
ENFS and SNFS stray light can be evaluated and subtracted point by point. 

All SALS instruments must include some device to stop the transmitted 
beam, like in ONFS. Moreover, the solid state sensors must be accurately posi- 
tioned with reference to the focus of transmitted beam. In an industrial instru- 
ment, the position of both the beam stop and the sensors must be electronically 
controlled, in order to correct the deformations due to mechanical stress and 
temperature dilatations. On the contrary, ENFS has an extremely simple and 
robust structure, which docs not require any adjustment. This makes ENFS 
suited for industrial applications of light scattering measurement, like particle 
sizmg. In Chapter [H we will show some particle size measurements performed 
with ENFS. 

In Chapter [51 we will show that SNFS can be used to make a measurement 
of the power spectrum of the non equilibrium fluctuations in a free diffusion 
process. Such fluctuations have never been observed by SALS, since they are 
extremely weak and involve mainly low wavevectors. This shows that NFS 
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Figure 3.7: Example of dark field image. 



teclmiques can be used to obtain measurements not possible with SALS. 

3.13 Why using NFS instead of classical microscopy? 

The near scattered field keeps only one feature associated to the observed sam- 
ple, the correlation function. For example, Figure [3771 shows a dark field image, 
and 13.81 shows the corresponding NFS image. The correlation function of the 
two fields is the same, but looking at the second immage we cannot figure that 
it comes from a set of discs. With near field scattcrig we will never distinguish 
an amoeba from a Paramecium: it is not a microscopy technique. Why using 
NFS instead of dark field? 

The first answer comes from the analysis of Figure 13.71 and 13.81 In Figure 
13.71 some of the discs arc in the focal plane, other aren't. If we want to analyze 
a dark field image, we must be able to select the particles which are in the focal 
plane, and exclude from the analysis all the others. On the contrary, NFS gives 
informations which are never affected by the misfocusing z: it provides three 
dimensional informations, and works well for thin samples as well as for thick 
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Figure 3.8: Example of Near Field Speckles ima 
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ones. 



If we want to analyze, for example, a colloid by a microscopy technique, 
we must use a thin sample, in order that the particles can be focused. If the 
concentration is low, it could be hard to find even one particle. Generally, one 
microscopic image could show only some particles. On the contrary, NFS can 
work on thick samples. We can use a given colloid, with any concentration, and 
put it in a cell so thick that it shows a suitable attenuation. 

Another reason leads to use ONFS technique instead of a dark field tech- 
nique. Dark field image intensity is given by equation 13.291 every calculation 
based on dark field images will concern the square value of the refraction index 
fluctuations. In facts, two point correlation functions of the images will rep- 
resent four point correlation functions of the refraction index fluctuations. We 
consider a fluid, for which the refraction index fluctuations has a give distribu- 
tion. We are interested in measuring the two point correlation, but dark field 
images allow us to work only on four point one. Of course, four point correlation 
function involves the two point one, but has a non trivial connected contribu- 
tion. In NFS images, every connected term in correlation functions vanishes: we 
can measure quantities directly connected with two point correlation functions. 

3.14 Meaning of the light path correlation func- 
tion. 

We have shown that NFS allows the measurement of the light scattered in a 
quite wide range of angles. If Reyleight Cans condition is met, the scattered 
intensity represents the power spectrum of the sample, evaluated in the trans- 
ferred wavevector. For the scattering at small angles, the spectrum is evaluated 
in the direction ortogonal to the incident beam. A measurement of this com- 
ponet of the spectrum leads, through a Fourier transform, to the correlation 
function of the light path through the sample: 



This quantity is directly accessible from NFS measurements. Its Fourier trans- 
form is the power spectrum for = 0: 



Through a measurement of the scattered light wc can know the power spec- 
trum in the plane perpendicular to the direction of the incident beam. This 




(3.75) 
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\5n{q,q,^Q)\' (3.76) 
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means that the Hght path correlation function bears less informations than the 
refraction index correlation function, but the light path correlation function is 
connected to the refraction index correlation function: 



If the sample is isotropic, its power spectrum depends only on the modulus of 
the wave vector. If we know the light path correlation function, wc can evaluate 
its Fourier transform, extend it to the three dimensions, and then, appling again 
the Fourier transform, wc obtain the refraction index correlation function. This 
operation is generally performed by using the well known Abel transform. 




Sn [x, z) 5n (x + Ax, z + Az) dxdzdAz 
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Chapter 4 

The experimental system. 



The whole system is sketched in Fig. 14.11 It consists in the light source, the 
sample cell and the image forming, capturing and processing system. In the first 
sections we describe all these parts and the criteria we followed to build our NFS 
instruments. In Sections 14.71 and 14. 81 we describe the three NFS instruments we 
buih. 

4.1 The optical system. 

The optical system is mainly built using elements supplied by Newport. The 
optical table is a VII3048 IsoStation. All the elements are mounted on X26 
rails, by using CN26 carriers. A picture of the system is shown in Fig. 14.21 

For the experiment descibed in Chapter [51 the optical axis must be vertical. 
The X26 rails are held in vertical position by mounting them on an X95 rail 
with suitable carriers: see Fig. 14.31 

4.2 The Ught source. 

The light source we use is a LED laser, coupled to a single mode fiber, supplied 
by Newport. In Tab. 14.11 we show the data supplied with the LED. Figure 
shows the LED laser and the fiber. 

The LED is housed in a Newport 700P temperature controller mount, con- 
nected to a temperature controller and a current driver. Since the output wave- 
lenght depends on the temperature, a temperature controller is needed. The 
LED laser has a built-in monitoring photodiode. The current driver monitors 
the output power to keep it constant. 

A single mode fiber is directly pigtailed to the LED. The other end of the 
fiber is terminated by an FC/PC connector. It is held by a Newport SL50BM, 
a gimbal mount for the adjustment of the azimuth and elevation, originally 
built for mounting mirrors or beam splitters. The beam output by the fiber is 
diverging; we coUimate it by sending it through a lens with a focal of 5cm. 
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Figure 4.1: The experimental system. 




Figure 4.2: Picture of the optical system. From left to right, we can see the 
adjustable mount that holds the fiber, the collimating lens, the cell, the micro- 
scope objective and the CCD camera. In the upper right corner there's the laser 
mount. 
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Figure 4.3: A view of the optical system for the experiment described in Chapter 
[9l All the optical elements are alligned in vertical direction. From the bottom, 
we see the optic fiber, held by an adjustable mount, the coUimating lens, the 
cell, held by the column on the left, the focusing lens, the blade, the neutral 
filter, and the CCD camera. 
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Item code 
Center wavelength 
Wavelength range 
Fiber output power 
Threshold current 
Operating current 
Rise/fall time 
Operating temperature 



LD-635-31A 
635nm 
zblOnm 
1.2mW 
10mA 
40mA 
1.5ns 

-10 to 40°C 



Fiber core/clad diameter 4/125//m - single mode 



Tabic 4.1: Data of the LED laser. 




Figure 4.4: The LED laser. 
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The resulting beam has a diameter of about 2cm. This diameter is enough for 
our experiments; in general, it must be selected as a function of the wavevector 
range [qmin, qmax] we want to measure. For ONFS and ENFS, D, the diameter 
over which the beam intensity is constant, must be selected in order to fulfill 
Eq. (13211) ■ The best choice is: 

D>50^ (4.1) 
a ■ 

For SNFS, the beam must have a diameter which covers a good statistical sam- 
ple: 

D > 20/qrmn- (4.2) 

Both the adjustment of the direction of the beam and its coUimation are 
not critical operations. The direction of the beam must be adjusted to hit a 
lens, centered at the optical axis, half a meter away from the fiber end. The 
coUimation is checked by measuring the beam diameter on a screen, near the 
lens and one meter away. The collimator is shown in Fig. 14.51 

The LED laser has been used in order to test its performance for industrial 
applications. The LED laser is much more compact and robust than a gas 
laser; it can operate immediately after it has been powered on and it does not 
generate too many heat. For these features it is ideally suited for industrial 
applications. Moreover, the overall cost of a LED laser device included in an 
industrial product can be made extremely low, as in the case of CD readers, 
though a laboratory LED laser equipement can cost as much as a classical gas 
laser. Moreover, the output of a gas laser must be spatially filtered before being 
used. A spatial filter is a critical component in an industrial equipement, since 
it must be extremely stable, and must be adjusted by micrometric actuators 
controlled by sensors, in order to correct the deformations due to heating and 
mechanical stresses. On the contrary, the single mode fiber output is more 
uniform than the output of a spatial filter, and requires no adjustment. 



4.3 The cell. 

The liquid samples we measure are held in a cell; the diameter D must be 
selected following Eq. gU) for ONFS and ENFS, or Eq. g21) for SNFS. 

For homogeneous samples, like colloids, the thickness must be selected in 
order to have a suitable attenuation of the main beam, about 1%. For ONFS 
measurements, the thickness of the cell and the volumetric particle density must 
fulfill Eq. (|3.63p : generally this condition is spontaneously met. 

For ONFS and ENFS measurements, the parallelism between the windows 
of the cell is not critical, nor the optical quality of them. Since the measured 
scattered light comes from different regions of the sample, we must provide that 
it is homogeneous. This implies that the thickness must be uniform, but an 
optical quality allignment is far beyond what is needed. On the contrary, SNFS 
requires optical quality windows: the well known "Foucault test" sees every 
deformation of the wavefront, no matter if the associated wavelength is long. 
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The cells we used are described in detail in Chapters [7] and [SI 



4.4 Objective, beam stop and blade. 

The objective must form an image of a given plane on the CCD sensor. The 
magnification M must be selected in order that the required wavevector range 
[(lmin,(lmax\ IS iuside the wavevector range the CCD sensor can measure: about 
[2 X 10^m"\2 X lO^m^i]. This means that: 



M ' M 



C [2 X 103m-i,2 X lO^m-i] (4.3) 



Moreover, the numerical aperture of the lens must be enough to resolve details as 
small as the smallest wavelength involved, 2'K/qmax, or, equivalently, to collect 
light scattered at an angle Qmax/k. 

In our experiments, we used a 20X microscope objective for high magnifi- 
cation factors, and an achromatic, 10cm focal length doublet for magnification 
factor around 1. An achromatic doublet has also been tested for high magnifica- 
tion factors, since we do not require the high quality of a microscope objective, 
nor an extremely wide numerical aperture. Experiments proved no different 
performances of the doublet compared with the microscope objective, but it 
was more difficult to obtain the required magnification. 

The objective lens must be placed so that it creates an image of a given plane 
on the CCD sensor. For ONFS and ENFS, the plane must be at a distance z 
from the sample fulfilling Eq. (|3.60p . The best choice is: 

z w 25^- (4.4) 

For SNFS: 

z < (4.5) 

^Qrnax 

For ONFS, the transmitted beam, focused by the objective, is stopped by an 
opaque or reflective element. In microscope objectives, the focal plane is inside, 
between two groups of lenses: we insert the beam stop through a hole. We tried 
three kinds of beam stops: a thin wire, a reflective wedge and an absorbing disc 
impressed by on a photographic film. The wire has a diameter of 70^m; it si 
stretched in the focal plane and is positioned by micrometric screws. It reflects 
the light inside the objective, and this could, in principle, increase the stray 
light. The photographic film we used are high contrast, black and white, 36mm 
photographic films. The beam stop is circular, but the beam is not completely 
blocked, thus increasing the stray light. The wedge was obtained by a steel 
blade; the edge was kept parallel to the optic axis. The upper part, in the 
direction from which the light comes, was cut at 45° and polished, in order to 
obtain a surface that reflects the main beam outside the lens mount, through 
a second hole. A section of the objective lens is shown in figure (|4.6p . This 
kind of beam stop is not symmetrical with respect to the optical axis. This 
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Figure 4.6: Section of the microscope objective and the beam stop. 

could increase the difficulty to process the data. During the experiments, all 
the methods showed to be almost equivalent. Figure [4771 shows the mount that 
holds the beam stop. 

For SNFS, a blade must be placed in the plane where the transmitted beam 
is focused. The blade must be extremely sharp: a razor blade is required. We 
mount it on a system with three micrometric screws, in order to accurately 
position it in the space. A picture of the Schlieren system is shown in Fig. 14.81 

4.5 CCD sensor. 

We used an industrial CCD camera: JAI CV M50. An image is shown in Fig. 
14.91 the data are provided in Tab. 14.21 

The output is a standard CCIR; since it is interlaced, an image is alwais 
formed by two fields acquired with a time delay of 0.04s, although the internal 
shutter allows to acquire a single frame in lOOns. 

The number of pixels and their dimension determine the wavevector range 
the CCD can directly measure: [2 x 10^m^^,2 x lO'^m^^]. Other wavevector 
ranges can be covered, by creating on image with a suitable magnification factor, 
but we cannot cover more than two decades. 

4.6 The acquisition and elaboration system. 

The frame grabber we used is an IC-RGB, from Imaging Technology. It performs 
an 8 bit digitalization of three standard composite video signals from the CCD 
cameras. It can be used to acquire simultaneously from two syncronized CCD 
cameras, for evaluating the intensity correlation function on two different planes, 
or to acquire from a single camera. 

The syncronyzcd acquisition from two cameras can be used to evaluate the 
three dimensional correlation function; the meaning of the three dimensional 
correlation function is explained in Appendix |Al 

The software for image elaboration was developed under Linux, written in 
C language. The drivers and the hbraries are the "IC-PCI" provided by "COM 
Optical Measurements Techniques" . The algorithms used to process the images 
are described in Chapt. [5]and[6l 



52 




Figure 4.7: Picture of the microscope objective with the beam stop. The beam 
stop is glued to the blue rod, held by a mount with three micromctric screws, 
for the adjustment of the position in every direction. 
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Figure 4.8: A view of the Schlieren system. From the bottom, we see the cell, 
the focusing lens, the blade, held by a micrometric mount, the neutral filter and 
the CCD camera. 
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Specifications 


CV-M50C 


Scanning system 


625 lines 25 frames/s 


CCD Sensor 


Monochrome 1/2" Hyper HAD IT CCD 


Sensing area 


6.6nimx4.8mm 


Picture elements effective 


752x582 


Elements in video out 


737x575 


Cell size 


8.6x8.3 


Resolution (horizontal) 


560 TV Unes 


Resolution (vertical) 


575 TV Unes 


Sensitivity 


0.5 lux, fl.4 


Sensitivity peak wave length 


500nm 


Wave length range 


400nm - 675nm (sensitivity > 50%) 


S/N ratio 


>56dB (AGC off, Gamma 1) 


Video output 


Composite VBS Signal 1.0 Vpp, 75 n 


Gamma 


0.45 - 1.0 


Gain 


Manual - Automatic 


Scanning 


2:1 interlace 


Accumulation 


Field - Frame 


Syncronization 


Internal, Xtal-generated - External 
HD/VD - random trigger 


HD Sync input - output 


4v, 75 n 


VD Sync input - output 


4v, 75 n 


Trigger input 


4v, 75 n 


Trigger input duration 


> HD interval 


WEN output (write enable) 


4v, 75 n 


EEN output (exposm-e enable) 


4v, 75 n 


Pixel clock out (optional) 


4v, 75 n 


Internal shutter 


Off,l/100s, 1/250S, 1/500S, 1/lOOOs, 
1/2000S, 1/4500S, 1.10000s 


Trigger shutter 


1/60,1/lOOs, 1/250S, 1/500S, 1/lOOOs, 
1/2000S, 1/4500S, 1.10000s 


Long time exposures 


one field to +oo. Duration between ex- 
ternal VD pulses. 


Operating temperature 


-5°C to +45°C 


Humidity 


20%-80% non-condensing 


Storage temperature 


-25°C to -t-60°C 


Storage humidity 


20%-90% 


Power 


12V DC ±10% 2.5W 


Lens mount 


C-mount 


Dimensions 


40mm X 50mm x 80mm 


Weight 


245g 



Table 4.2: Specification of Jai CV M50C camera. 
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Figure 4.9: Jai CV M50 camera. 




Figure 4.10: The optical setupt for the measurement of scattering from colloids 
some microns large with ENFS. 



4.7 ONFS and ENFS setup for colloid measure- 
ments. 



The systems are sketched in Figs. 14.101 and 14.111 

The measurements described in Chapters[7]and[8]cover the wavevector range 
hmn, Qmax] about [2 X 10'^m-\4 X lO'^m-i]. By using Eqs. (gHI) and (g^l), 
we obtain D ^ 5mm and z « 6mm. A 2cm beam diameter, obtained with lenses 
with 25mm diameter, as shown in Fig. 14.21 is enough to ensure that the intensity 
is constant over the length D. For the measurements described in Chapters [7] 
andlU we used larger lenses with 50mm diameter, with a larger beam diameter, 
in order to ensure a better uniformity, and z was increased acordingly. 

Following Eq. (|4.3p . we obtain the magnification: M ~ 20. We used a 20X 
microscope objective and numerical aperture of 0.45. 

The whole optical system for ENFS is shown in Fig. 14.21 

For ONFS, we insert a beam stop through a hole inside the lens mount. The 
beam stop and its adjustable mount is shown in Fig. 14.71 
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Figure 4.11: The optical setupt for the measurement of scattering from colloids 
some microns large with ONFS. 

4.8 SNFS setup for non equilibrium fluctuation 
measurements . 

The overall system is sketched in Fig. 14.121 

The range [qmin, qmax] of the fluctuations we measure is about [2 x 10^m~^, 2 x lO^m^^] , 
that is, the fluctuations range from ten microns to some millimeters. By using 
Eqs. ()4.2|) and (j4.5|) . we obtain D ^ 10mm and z < 125mm. The cell we used, 
described in Chapter [9l has an internal diameter of about 25mm. 

Following Eq. (|4.3p . we obtain the magnification: M = 1. Wc used an 
achromatic doublet with a 25mm diameter and focal length / ~ 100mm. To 
obtain the required magnification, p = q = 200mm. Since SNFS is affected 
by small inhomogencous fiuctuations of air temperature, we choose to put the 
collimating lens and the objective lens as close as possible to the cell, in order 
to prevent air movements. This resulted in a negative z. 

The whole optical system is shown in Fig. 14.31 



57 












Cell 



















Figure 4.12: The optical setupt for the measurement of non equilibrium fluctu- 
ations with SNFS. 
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Chapter 5 

ONFS data processing. 



In the chapter [3l we showed that, under some given conditions, the correlation 
function of ONFS images can be used to derive the intensity of the hght scattered 
by a sample. The calcultions were performed in the ideal case, in which the 
scattered light comes only from the sample. The presence of non ideal lenses 
and optical elements introduces an amount of undesired scattered light. This 
problem is common to every kind of scattering measurement; the undesired 
light, often referred to as stray light, is generally scattered at small angles. 

In standard scattering measurements, the effect of the undesired light is 
additive. It can be subtracted, since the stray light can be measured by a blank 
measurement. 

Dynamic scattering gives a way to distinguish the effect of the light scattered 
from elements that evolve with time from stationary ones. If the stray light 
comes from stationary elements, such as imperfections of the optical elements, 
its effect is to increase the correlation function, with no dependence on the 
delay. Thus the time dependent informations on the sample will be given by 
the bell shaped part of the correlation function, while the pedestal will contain 
informations on both the statically and dynamically scattered light. 

If a blank measurement is possible, a more refined subtraction of the stray 
light becomes possible [18], provided that the stray light constitutes a speckle 
field, that is, the field is gaussian. Such a data processing can be extended 
to ONFS too. In the following sections, we will find a way to subtract the 
effect of the stray light, first considering an unlimited number of images, taken 
at different times, and then a finite set of images. Then, we will describe the 
whole data processing algorithm. 

5.1 Effect of the stray light. 

From a set of ONFS images / {x) , we can measure the intensity correlation 
function: 

Ci {Ax) = {{I {x) I {x + Ax))} , (5.1) 
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where (•) is the mean over x, {■} is the mean over different images, and the 
intensity of the images / (x) is the intensity of the sum of SE {x), the field 
scattered by the sample, and Esl {x), the field of the stray light: 



So we obtain 
Ci (Ax) = 



Iix) = \dE ix) + EsL ix)\' 



\6Eix)f + \EsL {x)f + SE (x) E*sL {£) + 5E* (x) Esl (x) 



(5.2) 



\SE{x + Ax)f+ Esl (x + 



Since SE (x) is a random, circular gaussian field, the mean over different 
images of its odd powers vanishes; since Esl (x) is static, it can be considered 
as a costant, with respect to {•}, the average on the images: 

Ci (Af) ^ {(^\SE {x)f \SE{x + Ax)f)} + (^\Esl {xf \Esl {S + Ax)\^) + 
([\SE (f)|2} \EsL {x + Ax)f) + ^\EsL ix)\'{\SE {x + Ax)\^}) + 

{{6E (f) ,5^;* (f + Af )} E*sL (x) Esl{x + Ax)) + 

({SE* (x) SE (x + Ax)} Esl (x) E*^ {x + Ax)) . (5.4) 

The mean over the images {•} equals the mean over x, (•), for the field SE (x): 

Ci {Ax) = (^\SE {x)f \SE {x + Ax)\^) + (^\Esl {x)f \Esl {x + Ax)\^) + 
'\SE{x)\^)(\EsL {x + Ax)\^) + (\EsL ix)\') (\SE {x + Ax)'^' 



{6E (x) 6E* {x + Ax)) {E*sL (x) Esl{x + Ax)) + 

[SE* [x) 5E [x + A.f )) {Esl [x) E*sl [x + A.f)) . (5.5) 

Since both 5E (x) and Esl (x) are gaussian fields, we can use Siegert relation 
Eq. (j3.65p to express four-point correlation functions in terms of two-point ones. 

Ci{Ax) ^ (\SE{x)f'j (\6E {x + Ax)f'^ + \{SE{x)SE* {x + Ax))f + 
'^\EsL (f)l') {\EsL (f + A.f)|') + \{EsL {x)E*sL {x + Ax))\' + 
'\,SE (.f)|^) (\EsL {x + Ax)\') + (\EsL (x)l') (\SE {x + Ax)\') + 
{5E [x) 6E* {x + Ax)) {E*s^ {x) Esl [x + Af )) -I- 

{5E* {x) SE (x + Ax)) {Esl (x) E*sl (x + Ax)) . (5.6) 

We define {51) = (\5E {x)\^), {Isl) = (\Esl {S)\^), Cse (Ax) = {SE (x) SE* (x + Ax)), 



CsL (Ax) = {Esl (x) E*s^ {x + Ax)): 

Ci {Ax) = {SI)^ + \Cse {Ax)\^ + {Isl)^ + \Csl {Ax)\^ + 

2 {SI) {Isl) + Cse (Af) (Ax) + (Ax) Csl (Ax) . (5.7) 



60 



The result is that the stray light field correlation sums to the scattered field 
correlation: 

Ci (Ax) = ((«) + {IsL)f + \CsE (Af) + CsL (Af)l' . (5.8) 

In order to obtain informations about the correlation of the stray light field, 
we acquire a great number of images, with different scattered field, and we 
average them, thus obtaining the correlation function of the mean intensity 
{/ (a;)}. Then, we measure the correlation function of the mean intensity: 

q,} (Af ) = ({/ {x)] {I {x + Af )}) . (5.9) 

We evaluate the mean intensity {/(a?)} : 

{I{x)}^[\5E{x-)+EsL (x)r} = 

[\5E [x-t] + [\EsL {x)\^] + {5E (f) E*sL (f)} + {5E* (x) Esl {x)} . (5.10) 

Since Esl docs not depend on the image: 

{/ (f)} = {\5E [xf] + \EsL (f)l' + {5E {x)} E*sL (x) + {5E* (x)} Esl (x) . 

(5.11) 

Using the gaussian properties of the scattered light: 

{I{x)} = {SI) + \EsL{xf. (5.12) 
Now we can evluate the correlation function of the mean intensity: 

q,} (Af) = ([(,5/) + \EsL {x)\^] [{SI) + \EsL (f + Af)|']) = 
{Slf+{SI) (^\EsL {x;)\') + {SI) {\EsL (x + Af)|2) + (|i?sL (x)l' \Esl {x - 



Ax) I' 



(5.13) 



Using the gaussian properties of the field Esl ■ 

q,} (Ax) = {{SI) + {IsL)f + \Csl {Ax)f (5.14) 

From eq. (|5.12p . we can evaluate the mean value of the intensity of the 
images: 

{{I)} = {51) + {Isl) (5.15) 

Eq. (|5.8p . (|5.14[) . (|5.15[) give some informations about the field correlation 
of the scattered and stray light. If both the correlation functions are real and 
positive, the best evaluation of the field correlation function of the scattered 
field is: 

Ce (Ax) = ,JCiiAx)-{{I)}^ ^C[jyiAx)-{{I)f (5.16) 
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5.2 Correction for finite samples. 



In order to evaluate the correlation function of the mean intensity, we average 
a given amount of images, then we evaluate the correlation function of the 
obtained mean value. Since the number of images we average is finite, the 
correlation function will not correspond to that of eq. (|5.14p . For example, if 
the stray light vanishes, the mean intensity will still present fluctuations, due to 
the scattered light. These fluctuations vanish as the square root of the number 
of the averaged images, and consequently the correlation function becomes flat 
only for infinite samples. 

A similar problem arises when working with a stochastic, gaussian variable. 
If we have N values of the stocastic variable x, distributed with probability 
P (x) oc exp [— (x — xo) / (2cr^)] , we find that the best value for xq is the mean 
of the values x, and the best value for a is the root mean square displacement 
of the values x from xq. On the other hand, the average on a finite number of 
elements will be displaced from xq of an amount, vanishing as the square root 
of the number of the samples N, but so that the root mean square displacement 
of the data from the mean is alwais smaller than a. It is thus necessary to use 
the Bessel correction, dependent on the number of the samples N. 

Generally the Bessel correction is obtained in consequence of the "maximum 
likelihood" condition. This means that, given a set of values of a stochastic 
variable, and given a family of probability distributions, the parameters of the 
family must be selected in order to maximize the probability of finding the given 
data. Another approach is to find a suitable algorithm which gives the values 
of the parameters, from a set of data. The algorithm will be selected in order 
that the output values will be distributed around the true ones, with minimum 
square displacement. For a gaussian distribution, the two approaches give the 
same result. It is easy to show that, for example for a Heaviside distribution, 
the maximum likelihood condition fails to obtain the best results. 

In our case, the distribution function of the intensity is not gaussian. We 
will use weak condition, that is, we will look for an algorithm giving values 
which average to the true ones. In other words: we will try to avoid sistematic 
erroneous evaluations of the correlation function. 

We define ~ as the mean over N samples. In particular Cj is the correlation 
function of the averaged N images. To avoid sistematic errors, we must first 
evaluate {C/}: 



where the average intensity of N images / is given by the sum of the field 
scattered by the sample 6E (x) and the field of the stray light £^5 l (x) : 



{Cj{Ax)} = {{l{x)Hx + Ax))}, 



(5.17) 




(5.18) 



n 
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So we obtain: 



{Cj (Ax)} = — ^ i / {xf + \EsL (f)l' + {x) E*sl i'x) + SE^ (x) i?sL (x) 



iV2 

n.m ^ * 

(5.19) 

Wc can follow the calculations performed in Section [Ql to obtain Eq. ( [57 
In this case we obtain: 



|(5£;,„(.f + Ax)r 



EslIx 



{CriAx)} ^±Y.{\^En (f)l') {m-^ (f + Ax)r 

n.m 

l^J2\{6E,,{x)SE*^{x + Ax))\'+ 

[\EsL {x)f) {\EsL ix + Ax)f) + \{EsL {S)E*s^ {x + Ax))\^ + 
l^Y.{\^EAx)\'){\EsL {x + Ax)t 

n.m 

(\EsL (f)|')^^(|<5ii;„(f+Af)f 

n.rn 

^ {SEr. {x) SE*^ {x + Ax)) {E*sL (x) Esl {S + Ax)) + 

n.m 

^ ^ {SE: (x) SE,„, (x + Ax)) {Esl (x) (x + Ax)) . (5.20) 



The mean values can be calculated, provided that the two cases, n = m and 
n m are taken into account: 

{C'r (Ax)} = (61)' + 1 \CsE (Af)l' + (Isl)' + \Csl {Ax)\' + 2 (SI) (Isl) + 
^CsE (Ax) C*sL (Ax) + (Ax) Csl {Ax) . (5.21) 

The result reduces to eq. (|5.14p in the limit N ^ oo: 

{Cj (Ax)} = {{SI) + {Isl))'+^ \Csl (Ax)IV^ \Cse (Ax) + Csl {Ax)f . 

(5.22) 

From eq. (|5.8p . (|5.22p and (|5.15p we can evaluate the field correlation func- 
tion: 



CE{Ax) = JCi{Ax)-{{I)r- 



jJL_Cr{Ax) .^Cj (Ax) {(/)} 



2 

(5.23) 
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Figure 5.1: Near field intensity of the light scattered by a colloid of 5.2/im. 

5.3 Data processing algorithm. 

Once the experimental apparatus has been built, as described in Chapter [4l and 
the sample is placed in it, we acquire one hundred images. The electronic shutter 
of the CCD and its interlacement time must be so short that no evident evolution 
of the system happens during the exposure: for the samples we measured, that 
is colloids some microns large, with brownian movements, and non equilibrium 
fluctuations in the free diffusion of simple liquids, an interlacement delay of 
l/25s is sufficient. Moreover, different images must be completely uncorrelated. 
For a 10.0/im colloid, images must be grabbed at intervals longer than one 
minute, if only brownian movements are the source of decorrelation. For the 
non equilibrium fluctuations we studied, the interval was about one second. 

In figure [5TT] and [5?^ we show two typical images of the near field intensity of 
the fight scattered by colloids of 5.2/zm and 10.0//m. We can notice the different 
typical size of the speckles. 

For each image, we evaluate the correlation function. This operation is 
quite fast, since we can use a Fast Fourier Transform (FFT) algorithm. An 
FFT algorithm allows to evaluate the Fourier tranform of an M x TV matrix, 
with a number of arithmetic operations proportional to M N log (MN). By 
using Perceval relation, we can obtain the correlation function by making an 
FFT, evaluating the square modulus, and making an Inverse FFT (IFFT). This 
only requires a number of operation of the order of Af A^log {MN). By scanning 
every value of Ax, and averaging over every TV x M pixels, the number of 
operations would be of the order of (AfTV)^. Using FFT, care must be taken in 
order to correcly evaluate the correlations near the boundarys: FFT assumes 
periodic boundarys, so the image must be embedded in a bigger matrix, filled 
with zeroes. Since the FFT is faster if TV and M are powers of 2 ^9\, we used a 



64 



Figure 5.2: Near field intensity of tfie light scattered by a colloid of lO.Oyum. 



matrix of 1024 x 1024 points. After the correlation function has been evaluated, 
we normalize it, by dividing by the number of independent couples used to 
evaluate the correlation function. 

The correlation functions of every image are averaged, thus obtaining C/ (Ax). 
Fig. 15.31 and 15.41 show typical graphs of the intensity correlation function 
Ci (Aaf), for a colloid made of polystyrene spheres with diameters of 5.2^m 
and 10.0//m. We can notice that the correlation function has a maximum at 
Ax = 0, then decreases, until it reaches the plateau value, about one half the 
peak value. This behaviour is typical of every speckle field. 

Neglecting the stray ligth, we could evaluate the field correlation function 
by using the Siegert relation, Eq. (|3.65|) : 



where the mean intensity {(/)}, is obtained by averaging the measured intensity 
over every pixel of the image and over every image. In Fig. I5.5l and l5.6l are shown 
typical graphs of the field correlation function, calculated from the intensity 
correlation function, without any correction for the stray light. The correlation 
should vanish for Ax oo, in absence of stray light. 

In order to subtract the contribution of the stray light, we evaluate the 
correlation function of the average of all the images, thus obtaining Cj{Ax). 
The evaluation of the correlation function is obtained with the above described 
algorithm. In Fig. 15.71 and 15.81 are shown typical graphs of the correlation 
function of the mean intensity, for the two colloids. The graphs are not flat, due 
to the stray hght. 

Through Eq. (|5.23p we evaluate Ce (Ax), under the hypothesis that both the 
stray light field and the scattered light field have a real and positive correlation 




(5.24) 
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Figure 5.3: Intensity correlation function C/(A.t). for a colloid made of 
polystyrene spheres with diameter of 5.2/im 




Figure 5.4: Intensity correlation function C/(Aaf), for a colloid made of 
polystyrene spheres with diameter of 10.0/im 
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CE(Ax,Ay) (arbitrary units) 




Ax (m'^) 



x10'' 
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Figure 5.5: Field correlation function Ce (A^), not corrected for the stray light, 
for the colloid made of polystyrene spheres with diameter of 5.2/im 




Figure 5.6: Field correlation function Ce (Ax), not corrected for the stray light, 
for the colloid made of polystyrene spheres with diameter of 10.0/im 
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Figure 5.7: Correlation function of the mean intensity Cj (Ax), for the colloid 
made of polystyrene spheres with diameter of 5.2/im 
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Figure 5.8: Correlation function of the mean intensity Cj (Ax), for the colloid 
made of polystyrene spheres with diameter of 10.0/im 
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CE(Ax,Ay) (arbitrary units) 
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Figure 5.9: Field correlation function Ce (Ax) for the colloid made of 
polystyrene spheres with diameter of 5.2/im 



function. Typical field correlation function, corrected for the stray light using 
Eq. (|5.23p . are shown in figure and [OTH we can notice a signitificative 
increase in the smoothness of the graphs, with respect to Fig. 15.51 and 15.61 

We apply a Fourier tranform to the two dimensional correlation function 
Ce (Ax), thus obtaining the field power spectrum Se (q)- Since our samples 
are isotropic, we make an angular average of the power spectra, and represent 
our data as a function of the modulus q of q. The scattered intensity / (q) is 
then obtained by using Eq. (|3.14p . that is, simply relating each value of the 
power spectra, with wavelength g to a value of / (Q), where the relation Q (q) 
is given by Eq. ()3.13p . In Fig. 15. lll and 15.121 we show the measured I{q)- 



69 



CE(Ax,Ay) (arbitrary units) 




Figure 5.10: Field correlation function Ce (Aa?) for the colloid made of 
polystyrene spheres with diameter of 10.0/im 
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Figure 5.11: ONFS measurement of the scattered intensity of a 5.2/im colloid. 
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Figure 5.12: ONFS measurement of the scattered intensity of a 10.0/im colloid. 
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Chapter 6 

ENFS and SNFS data 
processing. 

In this chapter we describe the algorithm used to process ENFS and SNFS data, 
acquired by a CCD, in order to evaluate of the scattered intensity. From a set of 
images, taken at different times, we arc able to subtract the stray light, point by 
point: this is a noteworthy feature of the heterodyne techniques. The algorithm 
we describe is similar to the one which has been applied to shadowgraph, and 
should work for every heterodyne technique. 

The algorithm has been applied to SNFS, and the results are shown in Chap- 
ter [51 in this chapter, all examples refers to colloid measurements made with 
ENFS. 



6.1 Subtraction of the stray light 

As in every heterodyne technique, we measure the overall intensity / (x), gener- 
ated by the interference of an object field 6E (x) with a more intense reference 
beam with amplitude Eq. In our case, the object field is generated by the scat- 
tered beams, and the reference beam is the transmitted one. We measure the 
heterodyne signal i (x) : 

na;) = -J , (6-1) 

where /q is the mean intensity. At first order in 5E/Eq, Eq. (|3.22p holds: 

where we have assumed that i^o is real. Equation (|3.72p states that, under 
the conditions in which NFS works, the power spectrum of i (x) is the power 
spectrum of the electric field, the quantity we must measure in order to evaluate 
the scattered intensities. 
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We developed an algorithm to subtract the contribution of the stray light, 
directly on each image, point by point. This is a noteworthy feature of the 
heterodyne techniques, since in dynamic light scattering and in ONFS the sub- 
traction is possible only on the scattered intensity or on the correlation function, 
averages of square values. The scattered field can be decomposed into Es l (x) , 
the stray light field, and SE{x), the field of the light scattered by the sample. 
Both EsL {x) and SE (x) are much less intense than Eq, the reference field. At 
the first order, the resulting intensity is: 

/ (x) = E'^ + 2Eo^ [EsL {£)] + 2Eo^ [SE {x)] (6.3) 

In many cases, 6E (x) fluctuates in time and is correlated only on finite delays. 
On the contrary, stray light comes mainly from hard surfaces, and does not 
change as times go on. This is the case of the samples we studied. The spatial 
average of a scattered field is alwais wanishing; this property, along with the 
absence of correlation on different images, says that the average over many 
images of 6E (x) vanishes. 

In order to separate the contribution of the stray light, we average / (x) over 
many different images. Since the phase of SE (x) is random, its mean vanishes: 

{I{x)}^E^ + 2Eo^[Esl{x)]. (6.4) 

We use the symbol {•} for the mean over many samples, and the symbol (•) for 
the mean over x. The fluctuation I (x) — {I (x)} does not depend on Esl (x): 

I{x)^{I{x)}^2Eo^[SE{x)]. (6.5) 

Because of the conservation of the total intensity during the scattering process, 
by averaging {/ {x)} over the whole plane, we obtain /q: 

{{Iix)})=Ei=Io. (6.6) 

We can now evaluate the heterodyne signal i (x) , subtracting the the stray light 
contribution: 

{{Hx)}) 

6.2 Correction for finite samples. 

The quantity {I (x)} should ideally be evaluated by averaging infinite images. 
We obtain a good evaluation of it by averaging a great number N of images 
/„ (x), typically one hundred: 

From this evaluation, we obtain i (x) : 

^Axl = '-^^^^]j^ (6.9) 
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The average value I (x), evaluated over a given number of images, is sistem- 
atically different from the true mean value, in the direction that reduces the 
evaluation of the root mean square displacement from the mean. This problem 
is analogous to the one that leads to the so called Bessel correction for the eval- 
uation of the variance ct of a stochastic variable, from the knowledge of a finite 
number of stochastic values. 

We evaluate the correlation function of in (x) for each n, then we average 
them, thus obtaining Ci (Aa?). Now we want to evaluate {Ci (Ax)}, that is the 
mean value over infinite samples, in order to correct systematic errors: 

{Q (Ax)} _ ^ / /„ (^) - - 5] I,n (x) In (^+ Af) - - ^ /„. {a 

\I) y n=0 \ I m=Q J L m=0 

(6.10) 

The symbol (•) means the average over x. We can write/ (x) + SI {x) instead of 
In (x): 



{a (Ax)} 



TV ]v^( SIn{x)--Y.^^''-('^) <S/„(x + Af)--^«„(f+Af) U 

\^/ t n=0 \ L m=0 J L m=0 If) 



Evaluating the products: 



(6.11) 
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Since SI = 0: 
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n,m,/— 



iV- 1 1 



N 



{{SI {x)SI (x + A.f))} 



Now we can use Eq. (|6.5p : 



{C, (Ax)} 



N~ 1 
N 



Ce (Ax) 



{{5Im {x) Sin [x + Axmm) 



(6.15) 



(6.16) 



(6.14) 



The correlation function evaluated on N samples is proportional to the cor- 
relation function evaluated for N oo. The poportionality constant is the 
same of the well known Bcssel correction. 



6.3 ENFS data processing. 

Once the experimental apparatus has been built, as described in Chapter [4l in 
the absence of the sample, the CCD should be illuminated in a quite imiform 
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Figure 6.1: Background image. 



way. As a matter of fact, the illumination is never completely uniform, primarily 
because of the interference of the main beam with stray light. A typical image 
is shown in Fig. 16.11 We can easily see some sets of concentric circles, each due 
to reflections inside a lens, along with speckle patterns properly due to stray 
light. 

When the the sample is placed in the right position, we acquire about one 
hundred images for each measurement. The electronic shutter of the CCD and 
its interlacement time must be so short that no evident evolution of the system 
happens during the exposure: for the samples we studied, an interlacement 
delay of l/25s is sufficient. Moreover, different images must be completely 
uncorrelated. For a lOyum colloid, images must be grabbed at intervals longer 
that one minute, if only brownian movements are the source of decorrelation, 
while for the non equilibrium fluctuations we studied the images can be taken 
at intervals of Is. In figure 16.21 and 16.31 we show two typical ENFS images, 
generated by the interference of the main beam with the light scattered by 
colloids of 5.2/im and 10.0/im. The images show a mean intensity, modulated 
by the interference with the speckle pattern. We can notice the different typical 
size of the speckles. The set of concentric circles can be seen yet: the stray light 
will be removed with the following step. 

Once the images /„ (x) have been acquired, they are averaged, in order to 
evaluate I {x) and ^/). By using Eq. (|6.7p . we evaluate i{x), the heterodyne 
signal. Figures 16.41 and 16.51 show the heterodyne signal: since i {x) is negative, 
for some points, a constant intensity has been added. The images thus simply 
represent the ENFS images, cleaned from stray light and optical imprefections. 

The heterodyne signal of each image is then elaborated in order to obtain its 
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Figure 6.2: ENFS image of a 10.0/.(m colloid. 
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Figure 6.6: ENFS measurement of the field correlation function, for a 5.2/xm 
colloid. 



power spectrum. Simple Fourier transforming of the signal would be uncorrect, 
due to border effects. First of all, we evaluate the correlation function. This 
operation is quite fast, since we can use a Fast Fourier Transform (FFT) algo- 
rithm. An FFT algorithm allows to evaluate the Fourier tranform of an M x N 
matrix, with a number of arithmetic operations proportional to MiV log (MN). 
By using Perceval relation, we can obtain the correlation function by doing an 
FFT, evaluating the square modulus, and doing an Inverse FFT (IFFT). This 
only requires a number of operation of the order of MN log (MN). By scan- 
ning every value of Ax, and averaging over every N x M pixels, the number of 
operations would be of the order of (MN)'^. Using FFT, well known tricks can 
be used, in order to correct the boundary effects ^9\. Figure \6M and [6771 show 
the correlation functions thus evaluated. 

The correlation function evaluated following the above described algorithm 
suffers from shot and read noise, that is, for the noise due to the CCD light 
measurement and acquisition systems. Since such a noise is not correlated to 
the speckle field due to scattered light, the noise correlation function sums to the 
speckle correlation function. In order to evaluate the noise correlation function, 
we acquire a set of about one hundred images, before putting the sample in the 
system. Then, we apply the above described algorithm to the images, and obtain 
the correlation function of the noise signal. Figure 16.81 shows the correlation 
function of the noise signal. We can notice a marked peak in 0, quite narrow, 
representing the correlation inside a row, and a correlation between lines spaced 
by two pixels, due to interlacing. The correlation function of the noise signal is 
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Figure 6.7: ENFS measurement of the field correlation function, for a 10.0/im 
colloid. 
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Figure 6.8: Correlation function of the shot and read noise. 
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Figure 6.9: ENFS measurement of the scattered intensity of a lO.O/itm colloid. 



then subtracted by the overall correlation function. 

Once the correlation function has been evaluated, through an FFT we obtain 
the field power spectrum Se (q)- Since our samples are isotropic, we make an 
angular average of the power spectra, and represent our data as a function of 
the modulus q of q. The scattered intensity / (q) is then obtained by using Eq. 
(|3.14p . that is, simply relating each value of the power spectra, with wavelength 
g to a value of / (Q), where the relation Q (q) is given by Eq. ()3.13p . In Fig. 
I6.9l and l6.10l we show the measured / (g). 
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Figure 6.10: ENFS measurement of the scattered intensity of a 5.2/im colloid. 
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Chapter 7 



Performances of ONFS and 
ENFS on two colloidal 
samples. 



In this chapter we describe ONFS and ENFS measurements we made on two 
colloids and discuss the results. The optical setup has already been described 
in Sections [4.71 



7.1 The samples. 

The samples we used are monodisperse colloids made of polystyrene spheres 
suspended in water. In order to avoid sedimentation, we used a mixture of water 
and heavy water with a weigh fraction of about 0.5. The diameters of the two 
colloids we used are 5.2^m± 0.5^m and lO.O^mi 0.3^m, whose polydispersity 
is negligible. The diameters are quite large, since NFS gives advantages with 
respect to classical LS for small wavevectors. 

The colloids we measured are held in a cell with plane parallel windows. 
The diameter is about 4cm, since the sample and the beam intensity must be 
uniform on a length D, where D is given by Eq. p.64p . The thickness is about 
2mm. We selected the tickness and the particle concentration in order to have 
a suitable attenuation of the main beam, about 1%. For ONFS measurements, 
the thickness of the cell and the volumetric particle density are enough to fulfill 

Eq. (inisD. 

The liquid is held between the two windows by an 0-ring; the parallelism 
between the windows is not critical, nor the optical quality of them. Since the 
measured scattered light comes from different regions of the sample, we must 
provide that it is homogeneous. This implies that the thickness must be uniform, 
but an optical quality allignment is far beyond what is needed. 
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Figure 7.1: ONFS measurement of the 10.0/im colloid. A SALS measurement 
and a theoretical evaluation based on Mie theory are shown. 

7.2 Measurements. 

In Fig. 17. 11 17. 21 and 17. 31 are shown the results of ONFS measurement. The same 
samples have been analyzed also by a SALS instrument; the results are shown 
in the graphs. 

In Fig. 17. 4| 17.51 and 17.61 are shown the results of ENFS measurement. The 
same samples have been analyzed also by a SALS instrument; the result is shown 
in the graphs. 

Both ONFS and ENFS measurements agree with the SALS measurement, 
but ONFS gives less accuracy. 



7.3 What is the main source of error? 

The ONFS data processing is based on Siegert relation; the correlation function 
is evaluated by Eq. (|5.23p . First, we average the correlation functions of each 
image. The error in the evaluation of the intensity correlation decreases as the 
square root of the number of the samples, and this dependence is intrinsic in 
the stochastical nature of the technique. Then, we evaluate the field correla- 
tion function, by extracting the square root of the difference between the mean 
intensity correlation and the square mean intensity. Thus we obtain a quan- 
tity which converges to the field correlation function as the fourth root of the 
number of samples. 

On the contrary, ENFS gives directly the correlation function without any 
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Figure 7.2: ONFS measurement of the 5.2/xm colloid. A SALS measurement 
and a theoretical evaluation based on Mie theory are shown. 
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Figure 7.3: ONFS measurements of the two colloids. SALS measurements and 
theoretical evaluations based on Mie theory are shown. 
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Figure 7.4: ENFS measurement of the lO.O^m colloid. A SALS measurement 
and a theoretical evaluation based on Mie theory are shown. 



0.35 




0.0x10° 5.0x10^ 1.0x10^ 1.5x10® 2.0x10® 

q(m-') 

Figure 7.5: ENFS measurement of the 5.2/im colloid. A SALS measurement 
and a theoretical evaluation based on Mie theory are shown. 
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Figure 7.6: ENFS measurements of the two eoUoids. SALS measurements and 
theoretical evaluations based on Mie theory are shown. 

extraction of square root. We can notice that the plateau of Fig. I6.6l and r6.7l are 
more plain than in Fig. 15.91 and 15.101 The noise on the plateau is then spread 
across all the power spectrum when the Fourier transform is performed. 

The situation is worst than in dynamic light scattering. The reason is that 
we are working in two dimensions: in order that the power spectrum can be 
evaluated, the correlation function must approach faster than while in 

one dimension, it must be faster than We can perform an angular mean, 

but we gain only a term r~^/^. 

These considerations should explain why ONFS data are much less accurate 
than ENFS and SALS ones: the problem is the slow statistical convergence as 
the number of samples increases. In order to test this explanation, we performed 
some numerical simulations. We used the SALS data to simulate a power spec- 
trum. We created one hundred random fields, with gaussian probability and the 
given power spectrum. We obtained the homodyne and heterodyne signals, thus 
creating images similar to the ones we acquired during the experiments. Last, 
we processed the data with the above described algorithms. The simulations 
are only affected by the statistics: they are virtually free from any experimental 
error. In Fig. 17.81 and Fig. 17.71 we present, respectively, the results for the 
ONFS and ENFS simulations. They look quite similar to the corresponding 
experimental measurements of Fig. 17.61 and 17.31 This confirms that the main 
source of error, for our ONFS measurements, is the poor statistical quality of 
the samples. Since the quality increases as the fourth power of the number of 
the samples, we cannot make ONFS measurements better than ENFS ones, un- 
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Figure 7.7: Simulations of ONFS measurements of the two colloids. SALS 
measurements are also shown. 



less we process at least one million images: this is, for the moment, impossible. 



88 




Figure 7.8: Simulations of ENFS measurements of the two colloids. SALS 
measurements are also shown. 
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Chapter 8 

Particle sizing with ENFS. 



One of the most important applications of Light Scattering technique, from the 
industrial point of view, is particle sizing. Industrial particle sizers generally 
include some sensors, which measure the scattered intensity / (q), both at small 
and high angles. Generally, a mechanical system makes the powder or the colloid 
flow in a cell, so that a good statistical sample can be obtained. An algorithm, 
based on Mie theory, tries to find the distribution of particle diameters, which 
best fits the measured scattered intensity. 

In order to asses the reliability of ENFS applied to particle sizing, we an- 
alyzed some mixtures of two colloids. We prepared two colloidal solution of 
polystyrene spheres. In order that the density of the solvent matches the den- 
sity of the colloid, we used a solution of equal volumes of water and heavy water: 
the colloid was quite stable, and did not sediment evidently over some hours. 
The diameters of the two colloids are 5.2/im±0.5^m and 10.0/im± 0.3^m (sam- 
ples A and B). The refraction index of the polystyrene is 1.59, while the solvent 
has the refraction index of water, 1.33. Then, we prepared three mixtures of 
them, respectively with volume fractions of 1:1, 1:2, 2:1 of samples A and B. 
The scattered intensity was measured both with ENFS and a state-of-the-art 
SALS instrument. The data are presented in Figs. O [O [SJl KM KM 

We define a and /3 the volume fractions of samples A and B in each mixture; 
the scattered intensity of the mixture with a given a and (3 is la, [3 (q)- The 
scattered intensities la, [3 (q), obtained for the three mixtures, are compared 
with the scattered intensities Ia (q) and Is (q) of the two samples A and B. 
We evaluate the values of a' and /?' for which la^p (q) ~ a' I a {q) + P'Ib {q), by 
looking for the minima of the mean square deviation: 



The values of a' and /3' are the measured colloid concentrations, and must 
be compared with a and /3. Table [571] shows the measured values, a' and /?', 




E, Iliq) E, ^|(9)-[E, lA(q)lB{q)f 
E, Ia.,0{q)lB{q) E, /1(9)-E, Ia.,,3{q)lA{q) E, lA{q)lB{q) 

E, -f^-j) E, /|(9)-[E, lA{q)lB(q)Y 



(8.1) 
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Figure 8.1: Scattered light intensity measurement of a 5.2/im colloid (sample 
A). 
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Figure 8.2: Scattered light intensity measurement of a 10.0/xm colloid (sample 
B). 
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Figure 8.3: Scattered light intensity measurement of a mixture of the two sam- 
ples. Volume fractions: 1/2 A, 1/2 B 
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Figure 8.4: Scattered light intensity measurement of a mixture of the two sam- 
ples. Volume fractions: 1/3 A, 2/3 B 
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Figure 8.5: Scattered light intensity measurement of a mixture of the two sam- 
ples. Volume fractions: 2/3 A. 1/3 B 



Measured 


Actual 

a,(3 


0.52/1,0.465 
0.6M,0.32S 
0.31A,0.68B 


l/2A,l/2B 
2/3A, 1/3B 
1/3A,2/3B 



Table 8.1: Measured and actual values of volume concentrations of colloid A 
and B in the three mixtures. 



compared with the actual ones, a and [3. The agreement is quite good: this 
shows that ENFS is suited for particle sizing. 

The scattering data has been analyzed by an inversion algorithm based on 
Mie theory. Mie theory allows to evaluate the scattered intensity / (g) gener- 
ated by a given diameter distribution p (d) of dielectric spheres; the inversion 
algorithm looks for the distribution p (d) which gives the best approximation 
to the measured I{q)- The results are shown in Figs. 18.61 and 18.71 Two peaks 
are quite evident: they arc centered on the diameters of 5p,m and lOpm. The 
small peak centered around 8/im in the histogram of Sample A corresponds to 
the scattering of the dymers: the colloid is partially aggregated. The height of 
the peaks in Fig. 18.71 change accordingly to the fraction of the samples A and 
B in the mixture. 

It should be noticed that ENFS measures the intensity of the scattered beams 
with reference to the main beam. This allows to evaluate the particle concen- 
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Figure 8.6: Diameter distribution of the two colloidal samples measured by 
ENFS, obtained by an inversion algorithm based on Mie theory. The height of 
the bars is proportional to the intensity of light scattered by the particles in 
the range covered by the horizontal extension of the bar. Sample A is a 5.2fim 
colloid, and sample B is a 10/^m colloid. The two peaks are evident. Sample 
A shows a small peak centered around 8/im: it corresponds to the scattering of 
the dymers. 
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Figure 8.7: Diameter distribution of the mixtures of colloidal samples, measured 
by ENFS, obtained by an inversion algorithm based on Mie theory. The height 
of the bars is proportional to the intensity of light scattered by the particles in 
the range covered by the horizontal extension of the bar. The dotted curves are 
obtained by combining the values measured for samples A and B, shown in Fig. 
18.61 with coefRcients given by the volume fractions of the two samples. 
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tration, and not only the relative concentration of different particles. This is 
accomplished by using a single sensor; on the contrary, with SALS, the transmit- 
ted and the scattered beams must be measured by independent sensors, because 
the intensities are generally extremely different. This difference comes from the 
fact that SALS sensors measure the intensity of scattered beams, while ENFS 
measures the interference of them. For example, consider a sample that gener- 
ates a single scattered beam, whose intensity is 1/10000 than the transmitted 
one. For SALS, we need two sensors, one for measuring the scattered beam 
and one for the transmitted beam, and they require an accurate calibration. A 
single sensor could be used without calibration, but its dynamic range should 
cover 4 decades, and in this range it should be quite linear. For ENFS, the 
interference of the two beams generates a modulation of about 1/100. A single 
CCD array can easily measure such a modulation. 
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Chapter 9 



Non-equilibrium 
fluctuations in a free 
diffusion experiment 
studied with SNFS 

In this chapter, we describe a measurement of the power spectrum of the 
noncquihbrium fluctuations that arise during a free diffusion experiment. In 
Sect. 19.11 we discuss the origin of such fluctuations; in Sect. 19.21 we describe 
the physical system we studied; the results are shown in Sect. 19.31 The optical 
setup has been already described in Sect. 14.81 

9.1 Nonequilibrium fluctuations in free difl"usion 
processes. 

Diffusion is the fundamental mass transfer mechanism in many natural and 
technological processes. The diffusive transport can be interpreted by the sim- 
ple molecular random walk model. A more refined description requires the 
understanding of direct interaction between the diffusing particles and possibly 
hydrodynamic interactions. Both types of interactions may produce apprecia- 
ble changes in the magnitude of the effective diffusion coefficient D but, at any 
extent, diffusion is believed to give rise to an intimate and homogeneous remix- 
ing on matter. The general belief is that while the process occurs over quite 
microscopic distances, nothing peculiar should occur at any other lengthscale, 
except the molecular one where the random molecular diffusion takes place. 

It has been recently shown that, quite unexpectedly, gigant fluctuations are 
present during the diffusive remixing of two miscible phases of a binary mixture 
not too far from its critical point |13] . A fluctuating hydrodynamic description 
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has been developed [5D] , which indicates that gigant nonequihbrium fluctuations 
should be present during the diffusive remixing of fluids in general; moreover, it 
has been shown that the fluctuations can be considered the origin of the whole 
Fick flow m]. 

The presence of the fluctuations has been detected experimentally during the 
free diffusive remixing occurring in ordinary liquid mixtures and in macromolec- 
ular solutions [22[ 123] . The measurements concerned an ordinary, low molecular 
weight liquid mixture, an aqueous solution of a low molecular weight solid, a 
polymer solution and a protein solution, thus giving evidence that these anoma- 
lous fluctuations are a universal feature associated with spontaneous diffusion 
across a macroscopic gradient. 

A free diffusion experiment begins filling a cell with the two liquids, with 
the denser solution in the lower part to avoid convective instability. The two 
horizontal layers are initially separated by a fairly sharp meniscus. As soon 
as the two liquids came into contact, the diffusive remixing begins, and the 
meniscus rcpidly becomes smeared. The concentration profile inside the sample, 
initially a step function as a function of the height z, gradually evolves into an 
s-shaped function [24], until eventually, after a few days, the concentration 
becomes uniform throughout the sample. 

During the free diffusion process described above, intense fiuctuations arise. 
Their power spectrum S (q), with q in the horizontal direction, is given by: 

S{q) = So ^-—T. (9.1) 



The roll-off wave vector g^o is given by: 



where g is the gravity acceleration, ly is the kinematic viscosity, D is the diffusion 
coefficient, and f3 = ^ff^, T is the solutal expansion coefficient, a quantity that 
increases as increases the mismatching of the two liquids. The gradient Vc can 
be assumed roughly connstant in the region between the fluids, where diffusion 
takes place, and vanishes outside. The sample-dependent prefactor in Eq. (|9.1|) 
is given by: 

dn\^ Ac 



where Ac is the total concentration difference across the sample. 

The power spectrum S (q) dislpays a q^^ power low divergence at large wave 
vectors, q ^ qRO, ^md a saturation at a constant value at small wave vectors, 
9 ^ Qbo- The q^'^ power low is interpreted as the result of of the coupling 
of velocity fluctuations with concentration fluctuations, while the saturation is 
due to a stabilizing effect of gravity on long wavelength fluctuations. 

Moreover, the roll-off wave vector where the transition between the two 
regimes occurs gets smaller as /3gVc, and the the low wave vector value of the 
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power spectrum is roughly constant as free diffusion takes place, since the 
concentration near the upper and lower windows of the cell are initially constant. 

The nonequilibrium concentration fluctuations are originated from the cou- 
pling of velocity fluctuations with concentration fluctuations, due to the pres- 
ence of a macroscopic concentration gradient. This can be understood by simple 
naive arguments, discussed in detail in |20| and [25j . Suppose that a small parcel 
of fluid of linear size a undergoes a velocity fluctuation. This fluctuation will 
displace the parcel until the viscous drag will stop it in a time given approxi- 
mately by Tyisc = o?' /v^ V being the kinematic viscosity. If the displacement of 
the parcel occurs in a direction parallel to the macroscopic concentration gradi- 
ent, the parcel will be surrounded by fluid with different concentration. The life 
time of this concentration fluctuation is Tdifs = /D, and is much larger than 
the viscous time Tyisc^ as D v. Thus, in the presence of a macroscopic gra- 
dient, the effect of a short living velocity fluctuation is to induce a long lasting 
concentration fluctuation. Once a concentration fluctuation has been created, 
two mechanisms may contribute to its relaxation: diffusion and buoyancy. If 
the spatial extent of the fluctuation is small, then the fluctuation will soon dis- 
appear due to diffusion. This mechanism gives rise to the (?~* divergence of the 
static power spectrum at high wavevectors. As the wavevector increases, the 
velocity fluctuation lives for a shorter time, and can displace the parcel of a 
smaller amount, and this gives a factor q~^; moreover, the displaced parcel will 
be dissipated as q~'^. However, if the fluctuation is large enough, the buoyancy 
force acting on it will be able to restore the fluctuation in the layer of fluid 
having the same density in a time shorter than the diffusive one. This gives rise 
to the frustration of the g"^ divergence at smaller wavevectors. 

9.2 The cell. 

We have investigated the free diffusion process that takes place when two mis- 
cible fluids are brought in contact, the mixing between adjoining regions being 
kept as little as possible before a measurement sequence. The liquid sample we 
used was an aqueous solution of glycerol, with a weight fraction of 0.3. It was 
diffused into pure water. The cell was filled with the two liquids, with the denser 
solution in the lower part to avoid convective instability. The two horizontal 
layers are initially separated by a fairly sharp meniscus. As SNFS is an image 
forming technique, at least for big objects, we were able to thoroughly check 
the sample for spurious disturbances at the interface before starting collecting 
data. As soon as the two liquids came into contact, diffusion takes place, and 
the nonequilibrium fluctuations arise. 

The main difficulty is to fill the cell, keeping the interface between the two 
liquids as regular as possible. We used a Flowing Junction Cylindrical Cell 
(FJCC), a prototype developed for the study of nonequilibrium fluctuations 
in microgravity [26]. From Eq. (|9.2[) we see that the roll off depends on the 
intensity of _g, the gravitational acceleration. As g decreases, gravity acts at 
increasingly shorter wavevectors, and the divergence of fluctuations at small q 
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Figure 9.1: The Flowing Junction Cylindrical Cell, developed for measurements 
of non equilibrium fluctuations in free diffusion experiments in microgravity. 
The gray parts are made of perspex. The glass windows are blue in the drawing. 
The two liquids are injected through the holes: water in hole A and a solution 
of water and glycerol in hole C. They fill, respectively, the upper and the lower 
part of the cell. The two liquids come into contact in the middle of the cell, and 
the solution they form flows through the slit D and is extracted from the hole 
C. The porous rings, red in the drawing, make the flow more regular. 

becomes more evident. From Eq. ()9.3p we see that the intensity of the power 
spectrum for small wave vectors increases linearily in 1/g. This divergence of 
the intensity of fluctuations on g will be studied in an experiment performed 
on the Intarnational Space Station. A drawing of the prototype cell is shown in 
Fig. 19.11 a picture can be seen in Fig. 19.21 Two pipes feed the cell with the two 
liquids, in the present study water and a solution of water and glycerol, with 
a small pressure. The liquids enter in two ring-shaped chambers, from which 
they flow in the cylindrical cell passing throug porous elements. The flow is 
quite symmetric, due to the presence of the porous rings. The two liquids flU 
the cell, water on the top and glycerol on the bottom; they come into contact in 
the middle of the cell, and are pushed out the cell through a circular slit. The 
outgoing liquid is collected in a third chamber, passes through another porous 
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ring, and is collected by a pipe. 

The FJCC can be filled also in microgravity, since it is based on the flow of 
liquids. However, gravity greatly simplifies this task: since the denser fluid is in 
the lower part of the cell, big fiuctuations, created by macroscopic motions, relax 
due to buoyancy, while small fluctuations disappear quickly due to diffusion. 

9.3 Results. 

Although low angle light scattering techniques are very suitable to study long 
range correlated fluids, their sensitivity is hampered by the divergence of stry 
light at small wave vectors. The data presented in [521 US] were collected by 
using the shadowgraph projection technique. Shadowgraphy has traditionally 
been used to obtain a qualitative mapping of inhomogcncities in the index of 
refraction. However, very recently the technique has been reintroduced as a 
powerful quantitative tool to assess the features of long wavelength fluctuations 
in fluids [TSIITS]. 

The main problem of shadowgraph is the oscillatory behaviour of its transfer 
function: see Eq. (|3.37p . The scattering intensities around the zeroes of the 
transfer function cannot be measured: for example, in Fig. 3 of |22| . the values 
of S {q) are missing for q fa 2 ■ 10^m~^. Moreove, the region in which the 
zeroes are too frequent cannot be considered in the data analysis. The overall 
wavevector range covered about one decade. 

We measured the scatterig intensities at different times after the beginning of 
the diffusion process. The power spectra measured with SNFS are shown in Fig. 
19.31 They show the q~^ divergence and the saturation for small wavevectors. 
The roll off wavevector is about lO'^m^^, and is compatible with the value given 
byEq. (lO) 

The data cover about two decades in wavevectors: about ten times the range 
covered with shadowgraph. The quality must be compared with data shown in 
Fig. 2 of [13] , obtained with SALS on a similar, but quite peculiar system: the 
wavevector range covered by SNFS is slightly more wide. Moreover, it should 
be noticed that SALS gives no reliable results for the present system, that is 
for the nonequilibrium fluctuations in the free diffusion of simple fluids, due to 
stray light, since the scattering is weak and the wavelength associated to the 
process is quite long. The results we present are the best obtained for such a 
system up to now. 
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Figure 9.3: Measured power spectrum of non equilibrium fluctuations in the 
free diffusion process of glycerol in water. 
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Appendix A 



Three dimensional intensity 
correlation function. 



The speckle field generated by a stocastie sample is formed by speekles extending 
in both the orthogonal and parallel direction with respect to the direction of 
propagation of the wave. The intensity measured in a plane perpendicular to 
the direction of propagation varies as the plane is moved; small movements of 
the plane will give small variations in the intensities. As a matter of facts, the 
speckles appear and disappear as the plane is moved. This allows us to speack 
of the three dimensional appearence of the speckles. We will show that the 
speckles are elongated in the direction of the propagation of the light. If the 
diamer is a times a wave length A, their length is times A. 

In the following sections, we will show that the three dimensional correlation 
function of the intensity of the scattered light gives more informations than the 
two dimensional one; in some cases it is possible to determine the sign of the 
field correlation function, thus determining it completely. Moreover, in analogy 
to the quadratic relation between the diameter and length of a speckle, the 
longitudinal frequencies should be related to the square root of the frequencies of 
the sample: measuring the longitudinal correlations should double the dynamic 
of the system. 

A.l Evolution equation of the field correlation. 

For q ^ k, Eq. ()3.5|) can be approximated by: 



In this approximation, neglecting the phase term exjp{ikz), the field follows a 
Schroedinger equation: 




(A.l) 



(A.2) 
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The three dimensional field correlation is defined as follows: 

Ce ( Af , Az) = ^ / E{x,z)E{x + Ax,z^ Az) dxdz (A.3) 
^ J s 

In order to obtain an evolution equation for C (Aaf, Az), as Az increases, we 
evaluate the first derivative of the correlation function: 

If d 

Ce (Af , Az) ^ - I E (x, z) -^i—E {x + Ax, z + Az) dxdz (A.4) 
b J g oAz 



d 



dAz 
Using eq. (|A.2p 



9Az 



If i 

Ce {Ax, Az) = - / E{x,z) —\/'^E [x + Af, z + Az) dfdz (A.5) 
b Js 2fc 



The operator V acts on the first argument of E (x, z), thus it can be considered 
as acting on Ax: 

■^Ce [Ax, Az) = ^ V^-^ ^ E (x, z) E {x + Ax, z + Az) dxdz (A.6) 

This proves that the evolution equation for (Ax, Az), as Az increases, is a 
Schrocdinger equation: 

■ ^ Cb(Ax,Az) = -^V2C£(Ax,Az) (A.7) 



aAz ' ' 2k 

This equation can easily be solved in Fourier space: 

Ce {q, z) = Ce (g, z = 0) e"' 2fc (A.8) 
We can now extend eq. (|3.65p to the three dimensional case: 
C7(Ax,Az) = (/(x,z)/(x + Ax,z + Az)) = {if + \Ce {Ax, Az)\^ (A.9) 

A. 2 Gaussian speckles. 

In this section we consider gaussian speckles, and we evaluate their three di- 
mensional correlation function. 

Far field speckles are often generated by scattering a gaussian beam, so that 
the far field speckles have a gaussian correlation function. We consider gaussian 
speckles in near field, since the case is analitically solvable, and involves some 
calculations used in quantum mechanics. 

The field correlation function of the scattered light, in the plane orthogonal 
to z, is gaussian: 

Ax^ 

Ce {Ax, Az = 0) = Ce" 2a^ . (A.IO) 
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In the Fourier space: 



CEiq,Az = 0) = 2TT<7^Ce 2 



2 2 
cr q 



(A.ll) 



Using eq. (|A.8|) : 



_1^2 2 _ 



Coming back to real space: 



(A.12) 



Ce{x,z) = C 



(A.13) 



Now we evaluate the modulus of the field correlation function, the quantity 
needed in eq. (jA.9p to determine the intensity correlation function: 



9 ? 

X cr 



a^ + z'^/k'^ 

We can now evaluate the intensity correlation function for x = 0: 



and for z = 0: 



Ci (x, z = 0) = 



1 + ecr^ 



V 



(A.14) 



(A.15) 



(A.16) 



While the transverse correlation function follows a gaussian law, the longi- 
tudinal one is a Lorentzian, The diameter of the speckles is about a, while their 
length is cr^k. 



A. 3 Determination of the sign of the field cor- 
relation function. 

The power spectrum, that is Ce iq,z = 0), is real. If the sample is isotropic, 
it is symetric with respect to the origin, and then the correlation function 
Ce [x, 2 = 0) is real. The knowledge of the intensity correlation function with 
Az = gives the absolute value of the field correlation function. The sign of 
the field correlation function does not affect the intensity correlation function 
with Az = 0, but it can affect its value for Az ^ 0. 

In figure [ATT] and [A?2l we see an example of this effect. The figures show the 
graphs of the square correlation functions. The first is such that Ce {x, Az = 0) = 
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Figure A.l: Comparison between two square correlation functions. 

Ce {x, Az = 0) = sin {x) jx 



sin(a;) /x; in the second, the correlation function has the same absolute value, 
but alwais positive sign, for Az = 0. For Az = the square correlation func- 
tions are equal; their evolution for other values of Az are different. We can 
explain this fact considering the evolution of the positive and negative parts of 
the correlation function. The two parts evolve, and overlap, as Az increases. 
The interference of the two parts depends on the initial phase. 

The sign of the correlation function is alwais possible, in principle. The 
presence of errors could limit this possibility. 

A. 4 Longitudinal correlation. 

We want to derive the field correlation along the z axis. We consider its Fourier 
transform: 

Ce {Ax - 0, q,) = / Ce (q, z) e~'1-^AqAz (A.17) 

(27r) J 

Using eq. (jA.SP : 

.q^z 

Ce {Ax = 0, q,) = [ Ce (g, ^ = 0) e~* 2^ ~ *^^^dgdz (A.18) 

{2tt) J 
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Figure A. 2: Comparison between two square correlation functions. 
Ce (x, Az = 0) = |sin (x) /x\ 
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The integration over z gives a Dirac delta: 

Ce (Ax = 0, g,) ^ ^ Ce {q, z = 0) (5 (^|^ + iq,^ dq (A.19) 
In radial coordinates: 

Ce (Ax = 0, q,) ^ J Ce (g, 'P,z = 0) qS + iq,^ dqdip (A.20) 
If the sample is isotropic: 

Ce (Ax = 0, g,) ^ 2t: j Ce (<Z, ^ = (^|^ + dg (A.21) 
The integral can be evaluated: 

C£(Af =0,g,) = 27rC£;(\/2fc^,z = 0) v/2fc^ (A.22) 

The dynamic of an instrument measuring the longitudinal correlation func- 
tion is twice that obtained with the transversal one. This facts closely mirrors 
the quadratic relation between the diameter and the length of the speckles. 
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Appendix B 

Definitions 



Fourier transform: 



/(?) 



/ 



/ {x) e-'I'^dx 



(B.l) 



Inverse Fourier transform: 



(2^) 



1 



71 



I 



(B.2) 



Convolution: 

5(^) = 1/(^)1' 

NFS Near Field Scattering 

ONFS hOmodyne Near Field Scattering 

ENFS hEterodyne Near Field Scattering 

SNFS Schlieren-like Near Field Scattering 

LS Light Scattering 

SALS Small Angle Light Scattering 

IFS Intensity Fluctuation Spectroscopy 

CCD Charge Coupled Device 
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